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ABSTRACT 


The  defense  response  against  chemical  and  biological 
(Chem/bio)  weapons  has  gained  a  renewed  focus  in  light  of 
the  11  Sept  2001  terrorist  attack.  A  successful  response  to 
a  Chem/bio  attack  would  involve  measuring  and  predicting 
the  dispersion  of  a  toxic  cloud  in  the  atmosphere.  The  NPS 
Aeronautics  and  Astronautics  Department  is  working  together 
with  the  Meteorology  Department  on  a  technique  to  make 
toxic  cloud  measurements  using  an  Unmanned  Air  Vehicle 
(UAV)  .  In  support  of  this  mission,  the  UAV  will  require 
precise  and  accurate  air  data  (airspeed,  angle  of  attack 
["alpha"]/  and  sideslip  angle  ["beta"])  so  that  wind  data 
extraction  can  be  carried  out  from  air  and  inertial  data 
for  use  in  plume  dispersion  modeling.  The  efforts  in  this 
thesis  concentrate  on  the  air  data  system  to  produce 
precise  and  accurate  air  data  for  the  support  of  the 
Chem/bio  response  UAV  flights.  The  primary  concerns  are  the 
choice  and  design  of  the  air  data  system;  the  calibration 
of  the  system  using  the  flow  fields  from  computer 
simulation;  and  the  processing  of  air  data.  The  air  data 
extracted  will  be  used  for  wind  determination  so  that  the 
movement  of  the  Chem/bio  dispersed  agent  in  the  atmosphere 
can  be  predicted. 
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I  INTRODUCTION 


A.  BACKGROUND 

Chemical  and  biological  (Chem/Bio)  weapons  have  posed 
a  defense  response  security  concern  for  some  time  and  have, 
in  light  of  the  11  September  2001  terrorist  attack,  gained 
a  renewed  focus.  The  dispersion  of  a  toxic  cloud  is 
important  to  measure  and  predict  in  order  to  successfully 
respond  to  attacks  by  such  weapons.  The  work  in  this 

thesis  supports  the  Chem/Bio  response  project  which  is  a 
joint  project  between  the  NPS  Aeronautics  and  Astronautics 
Department  and  Meteorology  Department  to  demonstrate  and 
validate  a  method  for  the  synthesis  of  measurements  and 

predictions  to  aid  in  the  response  to  an  attack  by  chemical 
and  biological  weapons.  The  use  of  a  UAV  to  make  in  situ 

measurements  of  particulates  as  well  as  air  and  inertial 

data  for  winds  extraction  is  a  unique  contribution  to  the 
area  of  homeland  security. 

One  of  the  key  features  of  the  Chem/Bio  response 
project  will  be  the  movement  prediction  of  the  Chem/Bio 
plume  with  a  UAV  equipped  with  a  proper  sensor  suite.  In 
support  of  this  mission,  this  thesis  research  concentrates 
on  the  air  data  system  to  produce  precise  and  accurate  air 
data  for  the  support  of  the  Chem/bio  response  UAV  flights. 
The  primary  concerns  are  the  choice  and  design  of  the  air 
data  system,  the  calibration  of  the  air  data  system  using 
the  flow  field  from  computer  simulation  and  the  processing 
of  air  data  so  that  data  extraction  for  wind  determination 
can  be  carried  out. 
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B.  OBJECTIVES 

This  thesis  seeks  to  implement  an  air  data  system 
using  standard  sensors  on  the  NPS  FROG  UAV  (refer  to  next 
section  or  [4]  for  more  details  on  the  FROG  UAV)  .  With  the 
air  data  system  implemented,  test  flights  are  to  be  carried 
out  so  that  air  data  can  be  acquired,  processed  and 
analyzed.  Wind  data  can  then  be  extracted  for  support  of 
the  Chem/Bio  attack  response  UAV  project. 

In  addition,  this  thesis  explores  better  quality  and 
precision  air  data  sensors  for  use  on  the  FROG.  This  is  to 
improve  the  quality  of  air  data  so  as  to  enhance  the  wind 
data  extraction  which  is  highly  dependent  on  the  quality  of 
air  data  acquired. 

Details  of  both  of  these  objectives  will  be  covered  in 
the  subsequent  chapters . 
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II  AIR  DATA  SYSTEM 


A.  THE  FOG-R  UAV  (THE  FROG) 

The  NPS's  UAV,  the  FROG,  as  shown  in  Figure  II.  1,  has 
been  the  test  bed  for  advanced  control  and  airborne  sensor 
projects  at  NPS .  It  is  manufactured  by  BAI  Aerosystems  as 
the  BAI-TERN  (Tactically  Expendable  Remote  Navigator)  and 
derives  from  the  FOG-R  variant  of  the  BAI-TERN  used  by  the 
US  Army,  hence  the  name  'FROG' . 

It  is  a  small  high  wing  monoplane  with  conventional 
elevator,  rudder,  ailerons  and  flaps,  and  uses  servo-motors 
designed  for  radio-controlled  airplanes  to  drive  the 
control  surfaces.  It  has  a  pod  and  boom  fuselage  with  high 
mounted  wing,  sweptback  fin  and  rudder  and  low  set 
tailplane.  The  wings  are  fitted  with  flaps  and  the  engine 
is  mounted  above  the  wing  center  section.  The  landing  gear 
is  that  of  a  fixed  tricycle  type.  The  launch  and  recovery 
are  that  of  conventional  wheeled  take-off  and  landing 
respectively.  More  details  on  its  physical  characteristics 
and  engine  are  documented  in  Appendix  A. 
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Figure  II. 1  NPS  FROG  UAV 


Figure  I I. 2  FROG  UAV  3  View  Drawing 

B.  AIRSPEED  SENSOR 

Two  basic  pressures,  static  and  total  pressure,  are 
used  for  measurement  of  airspeed.  The  static  pressure  is 
the  atmospheric  pressure  at  the  flight  level  of  the 
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aircraft,  while  the  total  pressure  is  the  sum  of  the  static 
and  the  impact  pressure,  which  is  the  pressure  developed  by 
the  forward  speed  of  the  aircraft.  The  relation  of  the 
three  pressures  can  thus  be  expressed  by  the  following 
equation : 


pt  =  p  +  qc  (1) 

Where  pt  is  the  total  pressure,  p  is  the  static 
pressure,  and  qc  the  impact  pressure. 

In  incompressible  flow,  the  pressure  developed  by  the 
forward  motion  of  a  body  is  called  the  dynamic  pressure  q, 
which  is  related  to  the  true  airspeed  V  by  the  equation: 

q  =  \pV2  (2) 

Where  p  is  the  density  of  the  air  and  V  is  the  speed 
of  the  aircraft  relative  to  the  air. 

For  compressible  flow,  the  measured  impact  pressure  qc 
is  higher  than  the  dynamic  pressure  and  the  effects  of 
compressibility  must  be  taken  into  account.  Since  the  FROG 
operates  in  the  low  subsonic  range,  compressibility  effects 
are  ignored. 

The  airspeed  of  the  FROG  is  computed  based  on  dynamic 
pressure  measurement  using  a  pitot-static  probe  (mounted  at 
the  wingtip)  and  pressure  transducers.  The  dynamic  pressure 
is  'fed'  into  a  pressure  transducer  that  in  turn  converts 
it  to  analog  voltage  signal.  With  proper  calibration  and 
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application  of  Equation  (2),  the  airspeed  of  the  UAV  can  be 
computed.  Subsequent  sections  will  explain  the  rationale 
for  mounting  the  probe  at  the  wing  tip  rather  than  at  the 
nose  or  any  other  locations  on  the  FROG. 

The  pitot-static  probe  is  a  straight  conventional 
type.  The  total  length  of  the  probe  is  26  inches  long  and 
the  four  static  pressure  sensing  ports  are  located  1.125 
inches  aft  of  the  total  pressure  port;  6  diameters  (of  the 
probe)  aft.  Figure  II.  3  shows  the  pitot-static  probe 
mounted  on  the  aircraft  and  Figure  II.  4  shows  the  actual 
probe . 

The  pressure  transducer  used  is  a  0-4  inches  H2O 
differential  pressure  transducer.  It  takes  in  a  single 
power  supply  between  7.5  to  24  volts  and  gives  an  output 
signal  of  0  to  5  volts.  It  is  precision  temperature 
compensated,  calibrated  and  has  an  operating  temperature 
range  of  -10  °C  to  70  °C .  It  has  a  proof  pressure  of  10 
inches  H2O  and  hence  if  the  input  differential  pressure 
exceeds  the  proof  pressure,  the  sensing  element  of  the 
transducer  will  be  damaged.  Figure  II.  5  shows  the  actual 
pressure  transducer. 
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Figure  II. 3  Pitot-static  probe  mounted  on  FROG 


Figure  I I . 4 


Pitot-static  probe 


layout 
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Figure  II. 5  Pressure  transducer  for  airspeed  measurement 

C.  ANGLE  OF  ATTACK  AND  SIDESLIP  ANGLE  SENSOR 

Angle  of  attack  (alpha,  a  )  is  defined  as  the  angle 
between  the  relative  wind  in  the  plane  of  symmetry  and  the 
longitudinal  axis  of  the  aircraft.  Sideslip  (beta,  f3 )  angle 
is  defined  as  the  angle  between  the  wind  vector  and  the 
plane  of  symmetry.  Figure  II.  6  illustrate  the  definitions. 
Wind  vanes  mounted  on  potentiometers  (pots)  are  used  to 
measure  a  and  (3  .  Note  that  the  f3  vane  actually  measures 
'flank  angle  of  attack'  but  since  a  is  small,  true  [3  can 
be  approximated  by  'flank  angle  of  attack'. 
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Figure  II.  6  a  and  /?  definitions  (From  [8]) 

The  vane-pots  assembly  is  mounted  on  a  probe  that  is 
similar  to  the  pitot-static  probe  at  the  wing  tip.  The  wind 
vanes  are  attached  to  the  shaft  of  the  pots.  As  the  UAV 
pitches  and/or  yaws,  the  vanes  rotate  and  that  causes  the 
shaft  to  rotate.  The  rotation  of  the  shafts  changes  the 
resistance  of  the  pots  and  since  a  voltage  is  applied 
across  the  pots  from  the  FROG' s  power  supply,  an  analog 
output  voltage  signal  is  produced.  With  proper  calibration, 
the  a- ft  of  the  UAV  will  then  be  obtained. 

The  pot  used  is  a  low  torque  pot  that  has  a  maximum 
running  torque  of  0.05  ounce-inch.  It  has  a  resistance  of 
10  kQ  ±  5%  with  a  linearity  of  ±  1%.  The  maximum  power 
rating  is  2  watts  at  40  °C .  The  shaft  is  made  from  stainless 
steel  that  is  non-magnetic  and  non-reactive.  It  has  a 
rotational  life  of  500,000  shaft  revolutions  and  a  load 
life  of  900  hours.  It  is  to  be  operated  in  the  temperature 
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range  of  -55  °C  to  125  °C .  Figure  II.  7  shows  the  actual  a  -  /? 
vanes  mounted  on  the  pots.  The  data  specifications  for  the 
pots  can  be  found  in  Appendix  B. 

The  wind  vanes  that  will  be  mounted  on  the  pots  are 
made  of  light  weight  graphite  and  the  vanes  are  mass 
balanced  to  eliminate  biasing  errors. 


Figure  II.  7  a  -  (3  vanes  mounted  on  the  pots 

D.  DATA  ACQUISITION  PACKAGE 

Apart  from  the  sensors  described  in  the  previous 
sections,  the  remaining  air  data  system  components  include 
a  processor  with  12-bit  analog  to  digital  converter,  a 
spread-spectrum  radio  modem  (refer  to  Figure  II.  8)  and  a 
ground  station  computer  which  make  up  the  data  acquisition 
package.  The  processor  is  a  data  logger  operating  with 
either  TXBASIC  or  C  (programming  language) .  It  consists  of 
a  PIC  microprocessor,  a  +5  volts  switching  power  supply. 
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and  a  +3  volts  low  power  supply  for  energy  conservation. 
There  are  two  serial  RS232  ports  and  eight  TP  lines 
available  that  can  be  programmed  as  digital  serial  input, 
output  lines  or  other  functions.  The  communication  link 
between  the  UAV  and  the  ground  station  computer  is 
implemented  with  a  wireless  modem.  The  modem  uses  frequency 
hopping  spread  spectrum  technology  and  has  a  power  output 
of  1/3  watts.  It  is  capable  of  communicating  over  a  line  of 
sight  range  of  up  to  20  miles,  and  supports  data 
transmission  at  baud  rates  from  1200  bps  to  115.2  kbps.  The 
ground  station  computer  basically  does  the  real  time 
download  and  storage  of  the  air  data  transmitted  down  from 
the  FROG  via  the  modem. 


Figure  II. 8  Modems;  both  the  aircraft  and  ground  station 
computer  has  one  unit  each 
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E.  INTEGRATION  OF  AIR  DATA  SYSTEM 

The  schematic  of  the  air  data  system  is  as  shown  in 
Figure  II. 9.  For  clarity,  other  sensors,  for  example,  rate 
gyros,  humidity  sensors,  etc.  are  not  included  in  Figure 
II. 9.  The  air  data  system  is  such  that  the  data  acquired  is 
post  processed  rather  than  real  time  processed.  (It  can  be 
real  time  processed  and  real  time  processing  of  air  data  is 
researched  in  [9])  .  The  pitot-static  probe  allows  dynamic 
pressure  to  be  measured  via  pressure  tubing  by  the  pressure 
transducer.  The  transducer  basically  converts  the  pressure 
pulses  into  an  analog  signal  and  the  signal  is  transmitted 
to  the  processor.  As  for  a  -  (3  measurement,  the  pots 
generate  an  analog  signal  and  through  shielded  wires,  the 
signal  goes  into  the  processor.  The  processor  basically 
converts  both  analog  signals  to  digital  signals.  The 
digital  signals  then  go  into  the  modem  onboard  the  aircraft 


and  are 

transmitted 

to 

the 

ground 

station  computer. 

A 

MATLAB® 

script  file 

is 

used 

as  the 

software  package 

to 

acquire 

the  data  so 

that 

it 

can  be 

post  processed  either 

with  MATLAB®  or  any  spreadsheet  program. 


Figure  II. 9  Schematic  of  air  data  system 
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F.  PRECISION  OF  AIR  DATA  MEASUREMENT 

1 .  Airspeed  Precision  and  Accuracy 

a.  Total  Pressure  Measurement 

The  equation  of  airspeed  is  based  on  the 
measurement  of  the  impact  pressure  and  it  is  derived  from 
measured  values  of  total  and  static  pressures.  The  static 
pressure  at  successive  points  along  lines  of  airflow  past  a 
body  can  vary  widely,  whereas  the  total  pressure  along 

these  lines  of  flow  remains  constant.  For  this  reason,  the 
measurement  of  total  pressure  is  much  less  difficult  than 
the  measurement  of  static  pressure. 

When  aligned  with  the  flow,  almost  any  open  ended 
probe  registers  total  pressure  correctly  provided  that  the 
tube  is  located  away  from  any  boundary  layer,  wake, 
propeller  slip  stream,  or  engine  exhaust.  As  far  as  the 

FROG  is  concerned,  the  pitot-static  probe  is  mounted  at  the 
wingtip  which  has  minimal  effect  from  any  boundary  layer, 
wake,  propeller  slip  stream,  or  engine  exhaust. 

Jb.  Static  Pressure  Measurement 

The  flow  of  the  air  over  a  body  creates  a 

pressure  field  in  which  the  static  pressure  varies  from 
point  to  point,  while  the  total  pressure  at  all  points 

remains  the  same.  For  this  reason,  the  measurement  of 
static  pressure  is  much  more  complicated  and  therefore,  the 
problem  of  designing  a  static  pressure  measuring  system  is 
primarily  one  of  finding  a  location  where  the  static 
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pressure  error  varies  by  the  least  amount  throughout  the 
operating  range  of  the  aircraft. 

Figure  II. 10  from  [1]  shows  four  typical  types  of 
static  pressure  measuring  installations  (static  pressure 
tubes  ahead  of  the  fuselage  nose,  wing  tip,  vertical  fin 
and  fuselage  vent  on  the  side  of  the  fuselage) . 


p 

p 

p 

Ap 


free-stream  static  pressure 
local  static  pressure 

pressure  sensed  by  static-pressure  tube  or  fuselage  vent 
position  error,  p'-  p 


Figure  II.  10  Typical  installations  for  measurement  of 
static  pressure  on  an  aircraft  (From  [1]) 

(1)  Fuselage-Nose  Installation.  For  a 

given  position  of  static  pressure  orifices  ahead  of  the 
nose  of  a  fuselage,  the  static  pressure  error  depends  on 
the  shape  of  the  fuselage  nose  and  the  maximum  diameter  of 
the  fuselage.  Figure  II. 11  shows  the  variation  of  the 
static  pressure  error  with  the  blocking  effect  (reference 
to  Figure  II.  11,  x/D)  .  As  far  as  the  FROG  is  concerned,  if 
the  static  pressure  measurement  is  made  with  a  fuselage- 
nose  installation,  the  error  will  be  8%  as  the  blocking 
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effect  is  about  0.75  since  x  is  about  20  inches  and  D  is 
about  26  inches.  Using  a  longer  probe  does  reduces  the 
error;  however,  one  may  encounter  vibration  problems. 


D  =  Maximum  effective  fuselage  diameter 


Figure  II. 11  Static  pressure  error  variation  with  blocking 

effect  (From  ( 1 ] ) 

(2)  Wing-tip  installation.  For  a  given 

position  of  static  pressure  orifices  ahead  of  the  wing,  the 
magnitude  and  variation  of  the  error  depend  on  the  shape  of 
the  airfoil  section,  the  maximum  thickness  of  the  airfoil 
and  the  spanwise  location  of  the  boom.  The  magnitude  of  the 
errors  ahead  of  a  wing  tip  is  shown  in  Figure  11.12.  The 
figure  shows  that  the  error  is  highest  at  the  position 
closet  to  the  wing  and  it  decrease  rapidly  to  a  value  of 
about  1%  qc  at  an  orifice  location  of  x/t  =  10,  beyond  which 
further  reduction  in  the  error  is  minimal.  For  the  FROG, 
having  the  wing-tip  installation  means  an  error  of  about  1% 
since  x/t  =  8. 
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Figure  11.12  Static  pressure  error  at  various  positions 
ahead  of  wing  tips  (From  [1]) 

(3)  Vertical-fin  installation.  The 
factors  that  affect  the  measurement  of  the  static  pressure 
ahead  of  a  vertical  fin  are  similar  to  those  for  wing-tip 
installations.  As  far  as  the  FROG  is  concerned,  using  a 
vertical-fin  installation  to  measure  static  pressure 
presents  some  degree  of  mechanical  complexity  due  to  the 
aircraft  layout.  Hence  this  option  was  not  considered. 

(4)  Fuselage  vent  installation.  The 
fuselage  is  in  a  general  way  a  'static-pressure  tube' .  When 
the  fuselage  is  aligned  with  the  air  flow,  the  pressure  at 
a  vent  is  determined  by  its  location  along  the  body,  and 
when  the  fuselage  is  inclined  to  the  flow,  the  pressure  is 
dependent  on  the  radial  position  of  the  orifice  around  the 
body.  Similar  to  the  vertical  fin  installation,  the 
fuselage  vent  installation  is  not  considered  for  the  FOG-R 
due  to  the  complex  nature. 
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Considering  the  fuselage-nose  and  wing  tip 
installations,  the  latter  was  selected  due  to  the  very  low 
associated  errors. 

c.  Pressure  Lag  and  Leak 

The  pressure  at  a  sensor  can  be  different  from 
the  pressure  at  the  pressure  source  because  of  a  time  lag 
in  the  transmission  of  pressures.  The  pressure  at  the 
instrument  can  also  differ  from  that  at  the  pressure  source 
when  there  is  a  leak  in  the  pressure  system.  For  both 
cases,  the  sensor  indications  will  be  in  error  by  an  amount 
corresponding  to  the  pressure  drop  in  the  system. 

(1)  Pressure  lag.  When  the  pressure  at  the 
pressure  source  is  changing  rapidly,  as  in  the  case  of 
high-speed  dives  or  climbs,  the  pressure  at  the  instruments 
lags  behind  the  pressure  at  the  source  because  of  (1)  the 
time  for  the  pressure  change  to  propagate  along  the  tubing 
(acoustic  lag)  and  (2)  the  pressure  drop  associated  with 
the  flow  through  the  tubing  (pressure  lag) . 

From  [1],  errors  associated  with  acoustic 
lag  are  estimated  to  be  negligible  since  the  pressure 
tubing  used  in  the  FROG  is  relatively  short  (80  inches  of 
tubing  used)  .  As  far  as  pressure  lag  is  concerned, 
calculations  can  be  carried  out  to  find  the  lag  constant 
and  rate  of  pressure  change  and  pressure  drop.  However,  [1] 
concluded  that  for  relatively  simple  pressure  systems  with 
few  bends  and  tees  in  the  tubing,  which  the  FROG  certainly 
has,  pressure  lag  is  not  an  issue. 
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(2)  Pressure  leak.  The  pressure  at  the 

instrument  can  be  different  from  that  at  the  pressure 
source  if  there  is  a  leak  in  the  system  and  if  the  pressure 
outside  the  system  is  different  from  that  inside.  The 

magnitude  of  the  pressure  error  due  to  a  leak  depends  on 
the  size  of  the  leak  and  the  pressure  drop  across  the  leak. 

The  errors  in  airspeed  that  result  from  a 
leak  of  a  given  size  and  a  given  pressure  differential 
across  the  leak  can  be  determined  from  the  leak  rate  and 
lag  constant.  As  far  as  the  FROG  is  concerned,  it  is  valid 
to  assume  that  there  are  no  or  minimal  pressure  leaks  since 
very  little  pressure  tubing  connections  are  present  and 

numerous  checks  have  been  carried  out. 

2 .  a  and  /?  Precision  and  Accuracy 

Due  to  the  flow  field  created  by  the  aircraft,  the 

flow  angle  at  any  given  location  in  the  vicinity  of  the 
aircraft  will  generally  differ  from  the  true  a  and/or  /?  . 
As  the  flow  field  around  each  different  aircraft 
configuration  is  unique,  the  a  and/or  /?  errors  (i.e.  the 


difference  between 

local 

and 

true  angles) 

at 

a  given 

location  with  respect  to 

the 

aircraft  will 

vary 

from  one 

aircraft  to  another 

and 

more 

importantly. 

from 

one  trim 

flight  condition  to 

another . 

For  the  FROG, 

the 

position 

errors  will  be  estimated  via  numerical  computation  to  be 
presented  in  subsequent  chapters  and  sections.  The  three 
sensor  positions  which  have  been  proven  successful  and 
which  have  been  used  to  the  greatest  extent  are  positions 
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ahead  of  the  fuselage  nose,  ahead  of  the  wing  tip,  and  on 
the  forebody  of  the  fuselage. 

From  [2],  it  has  been  concluded  that  for  operation 

throughout  the  subsonic,  transonic  and  supersonic  speed 
range,  a  position  ahead  of  the  fuselage  will  provide  the 
best  installation.  For  the  FROG,  having  the  fuselage  nose 
installation  for  the  FROG  does  present  an  operational 

problem  since  the  nose  compartment  is  the  primary  payload 
compartment  with  many  electronic  components  in  it, 
requiring  frequent  access.  The  wing-tip  installation  was 
finally  chosen  due  to  ease  of  fabrication,  implementation 

and  operation.  For  a  and  /?  measurement,  the  wing  tip 

installation  does  present  a  greater  error  compared  to  the 
fuselage  nose  installation  due  to  the  lift  characteristic 
of  the  wing;  however,  that  can  be  corrected  (as  presented 
in  subsequent  sections) . 

G.  CALIBRATION  OF  SENSORS 

Calibration  of  the  sensors  is  necessary  to  correlate 
the  results  in  desired  units  with  the  output  digital  signal 
from  the  processor  in  the  data  acquisition  package.  The 
calibration  of  the  pressure  transducer  was  carried  out  and 
the  results  are  plotted  in  Figure  11.13.  For  this 
transducer,  the  output  voltage  was  found  to  vary  linearly 
with  the  delta  pressure  as  designed  by  the  manufacturer.  A 
first  order  curve  fitting  was  carried  out  and  yielded  the 
following  equation: 

Delta  pressure  [pascal]  =  0.0315  x  (output  voltage 

-  1.7781 


[volts ] ) 
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Figure  11.13  Calibration  curve  for  pressure  transducer 

Equation  (3)  allows  the  processing  of  the  airspeed 
data  acquired  to  be  carried  out.  The  pressure  range  of  the 
transducer  in  terms  of  inches  of  water  is  also  shown  in 
Figure  11.13  for  comparison  purposes.  Equivalent  indicated 
airspeed  will  be  obtained  from  Equation  (3)  and  this  will 
be  shown  in  Chapter  III. 

Similar  to  the  pressure  transducer,  the  a- (5  pots  were 
calibrated  and  the  results  of  the  calibration  of  both  of 
the  pots  are  presented  Figure  11.14. 
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Alpha  Beta  POT  Calibration  (Done  on  aircraft) 


Figure  11.14  Calibration  curve  for  a- ft  pots;  positive  /? 
is  defined  as  'wind  in  right  ear'’ 

From  Figure  11.14,  both  the  a  -  /?  pots  output  voltages 
were  found  to  vary  linearly  with  the  relative  angle.  The 
useful  range  of  the  pots  will  be  ±20°  since  the  FROG  is  not 
expected  to  fly  at  any  a  or  /?  greater  than  20°  or  lesser 
than  -20°.  Note  that  the  input  reference  voltage  supplied  to 
the  pots  was  4  volts.  A  first  order  curve  fitting  was 
carried  out  and  yielded  the  following  equations: 

a  [°]  =  -0.0099527  x  (Output  voltage  [Volts])  +  170.11 

(4) 

/?  [°]  =  0.0099514  x  (Output  voltage  [Volts])  -  129.2 

(5) 

Since  the  pressure  transducer  used  is  a  0-4  inches  H20 
differential  type  with  an  output  analog  signal  of  0-5 
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volts,  this  translated  to  a  full  range  of  0-40.4  meters  per 
second  of  equivalent  airspeed.  As  the  processor  in  the  data 
acquisition  package  has  a  12-bit  analog  to  digital 
converter  with  input  voltage  of  0-4.096  volts,  the 

40.4 

resolution  of  the  airspeed  sensor  is  — - - « 0.008  meters 

— 212 

4.096 

per  second.  Since  the  useful  range  of  the  a  -  (3  pots  (40°) 
does  not  span  an  output  analog  signal  of  0-5  volts  but 
rather  2.76-3.23  volts,  the  resolution  of  the  a  -  {5  sensor 

is  — -  ..  — ~~z~7 - «0.09°.  More  details  on  resolution  will  be 

3.23  2.76  12 

l4.096  4.096^ 

covered  in  Chapter  VI . 
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Ill  AIR  DATA  PROCESSING 


A.  FLIGHT  TEST  PROFILE  AND  PROCESSING  OF  AIR  DATA 

The  FROG  was  flown  at  Camp  Roberts  on  9  Oct  2002  to 
check  out  the  air  data  system  as  well  as  to  acquire  air 
data.  The  test  flight  was  conducted  at  the  McMillan 
Airfield  at  Camp  Roberts,  California.  The  flight  profile 
was  straight  and  level  passes  on  the  runway  heading  with  a 
turn  at  each  end  at  an  approximate  flight  level  of  50  feet 
AGL .  Figure  III.l  shows  the  flight  profile  of  the  test 
flight  obtained  with  the  GPS  unit  onboard  the  FROG  and  an 
aerial  photo  of  the  airfield.  The  flight  profile  is  in  GPS 
coordinates  (longitudinal  minutes  versus  lateral  minutes) . 


GPS  Lat  M  i  n 


Figure  III.l  Flight  profile  of  FROG  during  test  flight 
(left)  and  McMillan  Airfield  (right) 

The  air  data  acquired  from  the  test  flight  was  post 
processed.  Post  processing  of  the  air  data  involved 
filtering,  re-sampling,  calibrating  and  making  corrections. 
The  next  few  sections  explain  the  details  of  the  filtering, 
re-sampling,  calibrating  and  correction  process.  The  air 
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data  post  processed  will  be  analyzed  in  Chapter  IV;  the  air 
data  was  correlated  with  the  flight  profile  flown  to  check 
the  sanity  of  the  air  data. 

B.  FILTERING 

It  is  vital  to  filter  the  data  signals  to  remove  un¬ 
wanted  noise.  This  is  done  via  power  spectrum  density  (PSD) 
analysis  and  software  filtering  using  the  Signal  Processing 
Toolbox  (SPTool®)  in  MATLAB®.  The  airspeed  data  was  sampled 
at  a  frequency  of  20  hertz  while  the  a  - /3  data  was  sampled 
at  40  hertz.  Figures  III.  2  through  III. 4  show  the  PSD  plots 
of  the  three  signals.  It  is  obvious  from  the  plots  that  the 
air  data  of  interest  is  of  low  frequency;  note  the  peaks  of 
the  PSD  plots  at  very  low  frequencies.  The  frequency  of 
noise  was  in  the  range  of  1.5  hertz  and  above. 


Frequency  [hz] 


Figure  III.  2  PSD  plot  of  unprocessed  airspeed  data 
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Figure  III. 3  PSD  plot  of  unprocessed  a  data 


Figure  III.  4  PSD  plot  of  unprocessed  /?  data 

Using  the  MATLAB®  Signal  Processing  Toolbox,  low-pass 
filters  were  designed  (with  the  auto-design  feature)  using 
the  Chebyshev  type  II  HR  algorithm.  For  the  airspeed  data 
low-pass  filter,  the  cut-off  frequency  is  0.5  hertz  while 
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the  cut-off  frequency  for  the  a  -  ft  data  low-pass  filter  is 
2  hertz.  The  cut-off  frequencies  were  auto-generated  by  the 
auto-design  feature  and  manually  checked  for  sanity  before 
the  filters  were  applied.  Figures  III.  5  through  III. 10  show 
the  signals  before  and  after  filtering. 


Figure  III. 5  Unfiltered  airspeed  data 


Figure  III. 6  Filtered  airspeed  data  using  low  pass  filter 
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Figure  III.  7  Unfiltered  a  data 
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Figure  III.  8  Filtered  a  data  using  low  pass  filter 
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Figure  III.  9  Unfiltered  /?  data 
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Figure  III.  10  Filtered  (3  data  using  low  pass  filter 

Figure  III. 5  and/or  6  shows  the  general  idea  of  the 
test  flight  profile.  The  FROG  take-off  phase  starts  at 
about  near  the  5000th  data  point  and  six  runway  passes  with 
turns  were  made  (from  the  6000th  to  9500th  data  points)  .  More 
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details  of  the  flight  maneuvers  will  be  covered  in  Chapter 
IV. 

It  is  distinctive  from  Figures  III. 5  through  III.  10 
that  the  low  pass  filters  (with  the  specified  cut-off 
freguencies  as  mentioned  in  the  above  paragraph)  were  able 
to  reduce  the  level  of  noise  and  enhance  the  fidelity  of 
the  air  data. 

C .  RESAMPLING 

The  air  data  system  time  tagging  comes  from  the  GPS 
clock  which  is  streaming  data  at  a  frequency  of  10  hertz. 
In  order  to  ensure  that  analysis  of  air  data  is  possible, 
it  is  necessary  to  have  a  common  time  history  of  the  air 
data.  This  is  particularly  so  for  the  airspeed  and  a  -  (3 
data  since  these  data  were  sampled  at  different  frequencies 
from  the  GPS  clock  signal.  The  airspeed  and  a  -  f3  data  were 
resampled  using  the  MATLAB!i'  'interpft'  function  so  as  to 
achieve  a  common  time  history  of  the  data.  The  'interpft' 
function  is  a  one  dimensional  interpolation  using  the  FFT 
method.  Figures  III. 11  through  III. 13  show  the  time 
histories  of  the  filtered  airspeed  and  a  -  f3  data;  note  the 
common  time  base  of  the  data.  Figures  1 1 1. 5  through  1 1 1. 10 
show  the  x-axis  of  the  plot  to  be  data  points  whereas 
Figures  III. 11  through  III. 13  show  the  x-axis  of  the  plot 
to  be  time  now. 
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Figure  III.  11  Time  history  of  filtered  airspeed  data 
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Figure  III.  13  Time  history  of  filtered  /?  data 

D .  CALIBRATING 

As  mentioned  in  Chapter  II  section  G,  the  calibration 
results  (Equations  (3)  through  (5)  )  were  applied  to  the 
time  history  of  the  filtered  air  data.  Figures  III.  14 
through  III. 16  show  the  time  history  plots  of  the  air  data, 
filtered  and  calibrated  with  the  proper  desired  units.  Note 
that  the  airspeed  determined  is  equivalent  airspeed  since 
the  standard  sea  level  density  value  was  used  in  Equation 
(2)  . 
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Figure  III. 14  Time  history  of  filtered  and  calibrated 
equivalent  airspeed  data 


Figure  III.  15  Time  history  of  filtered  and  calibrated  a 

data 
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Figure  III. 16  Time  history  of  filtered  and  calibrated  /? 

data 

E.  STATIC  AND  DYNAMIC  CORRECTION 

The  airspeed  and  a  -  /?  data  have  to  be  corrected  for 
position  errors,  namely  AVpc,  Aapc  and  A  /?pc  (static 
correction)  and  for  induced  pitch,  roll  and  yaw  rates 
(dynamic  correction) . 

1.  Static  Correction  Through  Off-Body  Flow  Analysis 

The  position  error  corrections  were  carried  out  using 
numerical  simulation  with  a  software  suite.  Personal 
Simulation  Works  (PSW®)  from  [4],  which  is  comprised  of 
LOFTSMAN,  CMARC  and  POSTMARC.  PSW®  is  based  on  inviscid 
potential  flow  analysis  and  it  provides  for  panel  model 
development,  input  file  processing  and  the  visualization  of 
results.  Appendix  C  gives  details  of  the  software  suite 
PSW®.  Figure  III. 17  shows  the  numerical  model  of  the  FROG 

for  this  thesis  work.  This  numerical  model  was  generated 
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after  modification  of  a  previous  CMARC  simulation  carried 
out  by  the  author  of  [4]  so  that  position  errors  can  be 
computed  for  the  wing  tip  installations  of  the  pitot-static 
and  a  -  (3  probes.  The  CMARC  source  code  for  numerical 
simulation  can  be  found  in  [4] . 


Figure  III.  17  CMARC  numerical  model  of  the  FROG. 

CMARC  is  ideally  suited  for  off-body  flow  analysis. 
Note  the  off-body  streamlines  in  Figure  III. 17.  Off-body 

streamlines  may  be  placed  through  a  point  anywhere  in  the 

flow  field.  CMARC  will  then  follow  the  streamline  up  and 

downstream  the  distance  designated  in  the  input  file.  This 
is  particularly  useful  for  flow  visualization.  In  addition, 
CMARC  calculates  the  pressure  coefficient  and  velocity  at 
each  point  along  the  streamline.  For  this  thesis,  one 

streamline  was  placed  through  each  of  the  following 
locations:  static  pressure  source,  a  and  (3  vanes.  The 
pressure  coefficient  (from  CMARC)  is  used  to  calculate  the 
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static  source  position  error  while  the  velocity  vector 
(from  CMARC)  is  used  to  calculate  both  a  -  (3  probe  position 
errors  (a  strong  function  of  a  and  /?  respectively)  . 

a.  Analysis  of  Static  Source  Position  Error 

The  position  error  pressure  coefficient,  ACp-pc  or 
APp/qc,  is  a  function  of  free  stream  Mach  number  and  a 
provided  that  the  static  source  is  located  outside  of  a 
thick  boundary  layer  and  sideslip  is  minimized  (from  [5])  . 
In  the  case  of  the  FROG  which  is  in  incompressible  flow, 
APp/qc  becomes  only  a  strong  function  of  a  and  it  is 
basically  independent  of  airspeed  and  altitude. 

The  CMARC  numerical  simulation  was  carried  out; 
Figure  III. 18  shows  the  variation  of  the  position  error 
pressure  coefficient  against  indicated  a  .  The  second  order 
influence  of  a  is  clear  with  a  second  order  curve  fitting 
tightly  through  the  data  points. 


AOA  indicated  [degree] 


Figure  III.  18  Position  error  pressure  coefficient,  ACp-pc 
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Position  error  pressure  coefficient  can  be  turned 
into  position  corrections  for  airspeed  and  altitude.  The 
following  relations  were  developed  which  assumed  small 
errors  and  incompressible  flow: 


VfAC 

AF  =  — - p- 


pc 


and  AVpc  =VC-Vi 


(6) 


VAV 

AH  =- — al 


^ stdSo 


and  AHpc=Hc-Hi 


(7) 


Where : 

AHpc  =  altitude  position  correction 

AVpc  =  velocity  position  correction 

ACP  =  position  error  pressure  coefficient 

ostd  =  standard  day  density  ratio 

g0  =  gravitational  constant 

The  velocity  position  correction  was  applied  to 
the  airspeed  data.  Figure  III. 19  shows  the  variation  of  Vpc 
with  indicated  a  while  Figure  III. 20  shows  the  time  history 
of  filtered,  calibrated  and  statically  corrected  airspeed 
data  of  the  FROG.  Note  the  slight  difference  in  the 
magnitude  of  the  airspeed  in  this  figure  and  Figure  III.  14. 
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Figure  III. 19  Vpc  variation  with  indicated  a 


Figure  III. 20  Time  history  of  filtered,  calibrated  and 
statically  corrected  airspeed  data 

Jb.  Analysis  of  Alpha  Vane  Position  Error 

Local  flow  field  velocity  was  extracted  from  the 
off-body  streamline  listing  to  obtain  the  local  a  .  The  a 
vane  was  assumed  to  capture  the  x-z  component  of  the  local 
velocity  field  and  ignore  cross  flow  in  the  y  direction. 
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The  flow  field  velocity  was  turned  into  indicated  a  and  ex¬ 
position  correction  with  the  following  equations: 


a,  =  tan '(— )  and  Act  =a-a 

l  ^  ~V  *  pc  t  l 


(8) 


Figure  III.  21  shows  the  Aapc  as  a  function  of 
indicated  a  as  it  was  derived  from  CMARC  off-body  flow 
field  analysis.  The  a  correction  is  fairly  linear 
throughout  the  FROG  operating  envelope  and  the  corrections 
apply  at  all  incompressible  airspeeds  and  altitudes.  Figure 
III.  21  was  used  to  correct  the  a  data  of  the  FROG.  Figure 
III.  22  shows  the  time  history  of  filtered,  calibrated  and 
statically  corrected  a  data  of  the  FROG. 


AOA  indicated  [degree] 


Figure  III.  21  a  vane  position  error,  Aap_pc,  for  the  FROG. 
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Figure  III. 22  Time  history  of  filtered,  calibrated  and 
statically  corrected  a  data 

c.  Analysis  of  Beta  Vane  Position  Error 

Local  flow  field  velocity  was  similarly  extracted 
from  the  off-body  streamline  listing  to  obtain  local  f3 
angle.  The  beta  vane  was  assumed  to  capture  the  x-y 
component  of  the  local  velocity  field  and  ignore  vertical 
flow  in  the  z  direction.  Flow  field  velocity  was  turned 
into  indicated  /?  and  [3  position  correction  with  the 
following  equations: 

Pi  =  lan  '  (y)  and  A(3pc  =  J3t  -  (3 t  (  9 ) 

Figure  III.  23  shows  the  A  /3pc  as  a  function  of 
indicated  [3  .  The  (3  correction  is  fairly  constant 
throughout  the  FROG  operating  envelope  and  the  corrections 
apply  at  all  incompressible  airspeeds  and  altitudes.  Figure 
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III.  23  was  used  to  correct  the  /3  data  of  the  FROG.  Figure 
III. 24  shows  the  time  history  of  filtered,  calibrated  and 
statically  corrected  /?  data  of  the  FROG. 
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Figure  III.  23  /?  vane  position  error,  A  /?p-pc,  for  the  FROG 


Figure  III. 24  Time  history  of  filtered,  calibrated  and 
statically  corrected  (3  data 
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2. 


Dynamic  Correction  for  a  and  /? 


Since  the  a  -  (5  probe  is  not  located  at  the  CG  position 
of  the  FROG,  there  is  a  need  to  correct  the  a  -  (3  data  due 
to  induced  roll,  pitch  and  yaw  rates.  From  [6],  the 
correction  equations  are  as  follows: 


a 


—  K  ( r/  _ "a  n  i  -fa 

corrected  a  indicated  Tr  ^  T r 


V  V 


Z  n  A'  n 

*  corrected  ^  indicated  i  -  P  ^  ) 


V  V 


(10) 


Where : 

K a  and  K p  are  upwash  factors  on  a  and  /?  respectively 

p,  q  and  r  are  the  roll,  pitch  and  yaw  rates 
respectively 

x,  y  and  z  are  the  positions  of  the  a  -  /?  probe 
relative  to  the  aircraft  CG  in  accordance  with  the 
coordinate  system  as  shown  in  Fiqure  III. 25. 
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Figure  III.  25  Coordinate  system  for  dynamic  correction  to 

air  data  (From  [6]) 

The  application  of  the  upwash  factors  is  debatable  and 
the  effect  is  of  second  order  (as  mentioned  in  [6])  .  As 
such,  for  this  dynamic  correction,  a  value  of  1  was 
selected  for  Ker  and  K /?  .  With  the  dynamic  corrections  being 
carried  out,  the  values  of  a  - /?  are  those  of  true  values. 
Figures  III. 26  and  III. 27  show  the  time  histories  of  the 
complete  and  processed  true  values  of  a  and  /? 
respectively . 
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Figure  III. 26  Time  history  of  true  a 


Figure  III. 27  Time  history  of  true  /? 

3.  Density  Correction  for  Airspeed 

The  last  required  air  data  correction  is  the  density 
correction  for  airspeed  since  density  varies  with  location 

and  altitude.  The  meteorological  data  from  the  weather 

43 


balloon  showed  that  density  did  not  varied  much  from  ground 
level  to  the  altitude  of  the  flight  profile  on  the  day  of 
the  test  flight  (since  the  FROG  was  flying  at  low 
altitude) .  The  average  density  determined  from  the 
meteorological  data  was  1.182654  kilograms  per  cubic  meter 
and  this  value  was  used  to  correct  the  processed  airspeed 
to  get  the  final,  complete  true  airspeed.  Figure  III.  28 
shows  the  true  airspeed  of  the  FROG. 


Figure  III. 28  Time  history  of  true  airspeed 

F.  WIND  DATA  EXTRACTION 

As  mentioned  in  the  objectives,  the  motivation  is  to 
have  a  successful  response  to  a  Chem/bio  attack  which  in 
turn  involves  measuring  and  predicting  the  dispersion  of  a 
toxic  agent  in  the  atmosphere.  The  dispersion  of  the  toxic 
agent  is  dependent  on  the  wind  velocity  in  the  atmosphere. 
Hence,  wind  estimation  will  enable  dispersion  prediction  of 
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the  toxic  agent.  Wind  estimation  is  based  on  the  following 
fundamental  velocity  vector  equation: 


•v  =  lv  - 

w  v  b,i 


',RK 


(ll) 


Where : 


lVw  is  the  wind  vector  in  the  inertial  coordinate 
system 

lVbi  is  the  FROG  velocity  in  inertial  coordinate  system 
obtained  from  the  GPS  unit  onboard  the  FROG 


Ji?  is  the  rotation  transformation  matrix  from  wind  to 
inertial  coordinate  system 

Vbw  is  the  FROG  velocity  in  the  wind  coordinate  system 
obtained  from  the  pitot-static  probe 

A  Simulink®  model  (from  MATLAB®)  was  created  in  [9]  to 
extract  the  wind  data.  The  Simulink®  model  fundamentally  is 
based  on  the  above  mentioned  equation  and  it  takes  in  air 
data  and  computes  the  necessary  t'R ,  which  is  the  rotation 
matrix  to  the  inertial  coordinate  system,  so  that  wind 
extraction  can  be  carried  out.  Figure  III. 29  gives  a 
preview  of  the  schematic  of  the  Simulink®  model. 
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Figure  III. 29  Preview  of  Simulink®  model 

(From  [ 9 ] ) 


for  wind  extraction 
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IV  AIR  DATA  ANALYSIS 


A.  ANALYSIS  OF  AIR  DATA 

1 .  True  Airspeed 

Figure  IV.  1  shows  the  true  airspeed,  GPS  ground  speed 
and  altitude  (MSL)  level  of  the  FROG  during  the  test  flight 
as  mentioned  in  Chapter  III  Section  A.  Specific  legs  during 
the  test  flight  were  identified  and  are  as  indicated  in  the 
figure.  The  magnitude  of  the  wind  vector  can  be  estimated 
from  this  figure.  The  difference  between  the  true  airspeed 
and  GPS  ground  speed,  i.e.,  the  difference  between  the  blue 
and  pink  lines,  is  the  estimated  ballpark  wind  vector 
magnitude;  which  in  this  case  shows  that  the  wind  magnitude 
during  the  test  flight  was  no  more  than  10  meters  per 
second.  This  is  incorrect  from  a  strict  point  of  view, 
since  the  actual  wind  vector  magnitude  has  to  be  computed 
as  mentioned  in  Chapter  III  Section  F. 

The  FROG  started  its  take-off  run  at  about  time  =  250 
seconds.  Prior  to  that,  the  FROG  was  on  the  taxiway  being 
checked  out  for  its  functionality.  From  the  altitude  curve, 
it  was  evident  that  the  FROG  lifted  off  from  the  runway  at 
time  =  253.3  seconds.  The  aircraft  started  to  climb  and 
made  a  right  turn  followed  by  a  left  turn  to  line  up  for 
the  first  pass  of  the  runway.  A  total  of  six  runway  passes 
were  made.  The  duration,  average  speed  and  altitude  (MSL) 
level  for  each  pass  was  approximately  10  seconds,  38  meters 
per  second  and  300  meters  respectively.  Each  time  after  a 
pass,  the  FROG  turned  left  to  line  up  for  the  next  pass.  In 
the  first  part  of  the  turn,  the  FROG  was  climbing  from  an 
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altitude  (MSL)  of  300  to  320  meters  and  losing  its  airspeed 
from  38  to  30  meters  per  second  at  the  same  time.  As  the 
FROG  completed  each  turn  (second  part  of  turn)  to  line  up 
for  the  runway  pass,  it  descended  from  an  altitude  (MSL) 
320  to  300  meters  and  regained  its  airspeed  back  to  38 
meters  per  second . 
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Figure  IV. 1  True  airspeed  correlated  with  test  flight 

profile 

The  pilot  of  the  FROG  did  confirm  that  the  FROG 
was  indeed  performing  the  above  mentioned  maneuvers.  Hence, 
the  airspeed  data  was  concluded  to  be  sensible  and  valid. 
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2. 


True  a-p 


Figure  IV. 2  shows  the  true  a-p  of  the  FROG.  Note  that 
the  figure  only  shows  the  true  a~P  after  the  FROG  had 
lifted  off  from  the  runway,  i.e.,  from  time  =  250  seconds, 
for  the  sake  of  clarity.  The  runway  pass  is  as  indicated  in 
the  figure.  The  first  observation  to  be  made  from  the 
figure  is  that  the  a  was  always  positive  throughout  the 
time  which  the  FROG  was  airborne.  This  makes  sense  since 
the  FROG  was  not  expected  to  fly  at  negative  a  .  The  second 
observation  is  that  the  P  fluctuated  between  ±  5°  most  of 
the  time.  This  also  makes  sense  as  the  pilot  was  constantly 
making  rudder  and/or  aileron  input  to  prevent  the  FROG  from 
flying  in  a  'crab  like'  manner. 


Figure  IV. 2  True  a-p  correlated  with  test  flight  profile 
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As  mentioned  in  the  previous  section,  the  FROG  turned 
left  in  between  runway  passes  and  each  turn  was  comprised 
of  a  climb  with  decreased  airspeed  followed  by  a  descent 
with  increased  airspeed.  As  the  FROG  initiated  each  left 
turn  (climb  and  decreased  airspeed) ,  the  a  decreased  from 
approximately  10°  to  5°  while  the  /?  increased  from 
approximately  0°  to  5°.  The  completion  of  the  left  turn 
(descent  and  increased  airspeed)  was  achieved  with  a 
increasing  from  5°  to  10°  and  /?  decreasing  from  5°  back  to 
0°.  Since  /?  was  defined  as  positive  when  'wind  in  right 
ear'  which  implied  a  left  turn,  this  is  consistent  with  the 
data  presented.  Note  that  during  the  duration  of  the  runway 
passes,  [3  fluctuated  between  ±  3°  which  was  evidence  of  the 
FROG  compensating  to  fly  straight  along  the  runway  passes. 


Figure  IV.  3  Sample  close-up  of  true  a  -  ft 

Figure  IV.  3  shows  a  sample  close-up  of  true  a  -  /?  with 
the  data  points  from  time,  t  =  300  to  305  seconds.  Note  the 
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closely  'packed'  data  points,  indicating  that  the  data 
points  represent  actual  data,  not  high  frequency  noise. 
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V  POSSIBLE  IMPROVEMENTS  TO  AIR  DATA  SYSTEM 


A.  PRESSURE  TRANSDUCER  COMPARISON 

The  pressure  transducer  that  was  used  on  the  FROG  (for 
ease  of  identification,  it  is  named  as  transducer  A)  was 
compared  to  another  type  of  transducer;  transducer  B. 
Appendix  D  documents  the  detailed  data  specifications  of 
transducer  B.  A  wind  tunnel  experiment  was  conducted  to 
compare  both  transducers . 

1 .  Wind  Tunnel  Experiment 

The  wind  tunnel  experiment  was  conducted  with  a  pitot- 
static  probe  mounted  in  the  wind  tunnel  to  measure  dynamic 
pressure.  The  dynamic  pressure  was  given  as  input  to  both 
transducers  at  the  same  time  so  that  the  quality  of  the 
output  signals  from  both  transducers  can  be  compared.  Both 
of  the  transducers  output  signals  were  processed  by  a 
processor  with  a  12-bit  analog  to  digital  converter  which 
was  sampling  data  at  a  frequency  of  40  hertz.  The  first 
part  of  the  experiment  involved  operating  the  wind  tunnel 
at  different  speed  settings  while  the  second  part  involved 
operating  the  wind  tunnel  at  a  particular  speed  setting  but 
subjecting  both  transducers  to  vibration. 

2 .  Wind  Tunnel  Experiment  Results 

Figure  V.I  shows  the  raw  output  signal  of  both 
pressure  transducers  from  the  first  part  of  the  experiment. 
The  step-like  shape  of  the  curve  is  due  to  the  different 
wind  tunnel  speed  settings  carried  out  during  the  first 

part  of  the  experiment.  The  magnitude  of  the  signal  is  not 
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important  but  rather  the  fluctuation  of  the  signal  at  a 
particular  voltage  is  of  importance  since  that  translates 
to  the  noise  level  of  the  transducer  which  in  turn 
determines  the  quality  of  the  transducer.  Although  both 
signals  from  Figure  V.l  look  almost  the  same  in  terms  of 
shape,  upon  careful  scrutiny,  the  blue  signal  shows 
obviously  that  transducer  B  is  of  better  quality  as  the 
noise  level  or  signal  fluctuation  is  less  than  that  of 
transducer  A.  This  is  evident  from  Figure  V.2  which  shows  a 
close-up  of  the  output  signals  from  both  transducers  from 
the  50th  to  100th  data  point  at  the  same  scale.  A  PSD  plot  is 
carried  out  for  both  transducers  to  further  verify  the 
noise  level  and  the  frequency  content  of  both  signals. 
Figures  V.3  shows  the  PSD  plots  of  both  transducers. 


Figure  V.l  Output  signals  of  transducer  A  and  B  during 
first  part  of  wind  tunnel  experiment 
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Figure  V.2  Sample  close-up  of  output  signals  of  transducer 
A  (left)  and  B  (right)  from  50th  to  100th  data  point 

Figure  V.3  was  produced  from  the  first  part  of  the 
wind  tunnel  experiment  (different  wind  tunnel  speed 
settings) .  Note  that  Figure  V.3  shows  a  lot  of  peaks  in  the 
range  of  6  hertz  and  above,  evidence  of  the  fact  that  the 
noise  associated  with  the  transducers  are  of  frequencies  of 
6  hertz  and  above.  The  peaks  in  the  region  of  less  than  2 
hertz  for  both  plots  show  the  fundamental  frequency  of  the 
signal  (dynamic  pressure)  being  measured.  The  fact  that  the 
PSD  curve  for  transducer  A  is  of  a  higher  value  than  that 
of  transducer  B  points  to  the  fact  that  transducer  A  output 
signal  contains  more  noise. 
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Figure  V.3  PSD  plot  of  transducer  A  and  B 

The  second  part  of  the  wind  tunnel  experiment  involved 
analyzing  the  signals  from  the  transducers  by  subjecting 
the  transducers  to  external  vibration.  The  mass  elements 
inside  the  pressure  transducers  respond  to  vibration  as 
accelerometers,  another  noise  source.  The  vibrations  were 
generated  from  a  Siemens  S40  cell  phone.  This  was  done  to 
simulate  the  vibrations  that  the  transducers  would 
encounter  when  installed  onto  the  FROG.  Figure  V.4  shows 
both  the  output  signals  of  the  transducers  with  both 
transducers  subjected  to  vibrations  while  Figure  V.5  shows 
a  sample  close-up  of  the  signals  from  the  450th  to  the  500th 
data  point.  Note  the  obvious  noisy  signal  and  the 
relatively  large  signal  fluctuation  produced  by  transducer 
A. 


56 


Figure  V.4  Output  signal  of  transducer  A  and  B  during 
second  part  of  wind  tunnel  experiment 


Data  points  Data  points 


Figure  V.5  Sample  close-up  of  output  signals  of  transducer 
A  (left)  and  B  (right)  from  450th  to  500th  data  point 

Similarly  to  the  first  part  of  the  wind  tunnel 
experiment,  a  PSD  plot  was  carried  out  for  both  transducers 
to  further  verify  the  noise  level  and  the  frequency  content 
of  both  signals.  Figure  V.6  shows  the  PSD  plots  of  signals 
from  both  transducers  subjected  to  vibrations.  Comparing 
with  the  PSD  plot  in  Figure  V.3,  there  are  a  lot  more  peaks 
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this  time  due  to  the  fundamental  vibration  frequencies. 
Note  how  noisy  the  transducer  A  is  compared  to  transducer 
B.  Compare  the  peaks  from  the  range  of  more  than  2  hertz. 


Figure  V.6  PSD  plot  of  transducer  A  and  B 

Figures  V.l  through  V.6  clearly  suggest  the  fact  that 
transducer  B  is  of  a  better  quality  in  terms  of  noise 
production . 

B.  MULTI -HOLE  PROBE 

Fluid  mechanics  have  employed  a  variety  of  single  and 
multi-hole  probes  in  research  and  design.  The  most  common 
is  the  well-known  pitot  static  tube  which  is  used  on  the 
FROG,  while  another  single-hole  probe  is  the  Preston  tube, 
which  is  practically  a  f lattened-down  pitot  tube.  Multi¬ 
hole  probes  have  also  been  used  extensively  to  provide 
velocity  magnitude  and  direction  with  high  accuracy  for 
angles  as  high  as  75°  for  standard  five  and  seven-hole 
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probes,  and  165°  for  omni-probes.  Multi-hole  probes  are 
designed  to  determine  velocity  through  direct  measurement 
of  the  pressures  at  the  probe  tip  and  then  using  the 
pressures  to  calculate  the  velocity.  These  probes  provide 
point  measurements  of  the  flow  velocity  vector  and  the 
total  and  static  pressure  with  minimal  interference  to  the 
flow . 


Alternative  tools  for  point  measurements  in  a 
flowfield  are  the  hot-wire  anemometer  and  the  Laser-Doppler 
Velocimeter.  A  quick  comparison  of  the  advantages  and 
disadvantages  of  the  three  main  tools  are  as  follows 
(extracted  from  [7]): 

Hot-Wire  Anemometry 

A  hot-wire  anemometer  imposes  low  interference  on  the 
flow  and  provides  the  highest  possible  frequency  response. 
On  the  other  hand,  it  requires  often  and  tedious 
calibration  and  frequent  repairs  of  damaged  wires. 
Moreover,  to  generate  more  than  one  component  of  the 
velocity,  one  has  to  use  multiple  wires,  which  are 
sensitive,  fragile  and  hard  to  calibrate. 

Laser-Doppler  Anemometry 

This  method  relies  on  the  processing  of  laser  light 
scattered  by  micro-scale  particles  moving  with  the  flow. 
Since  the  technique  is  optical,  it  presents  no  interference 
with  the  flow.  It  does  not  require  calibration  but  the  user 
must  seed  the  flow.  Seeding  is  a  tedious  and  not  always 
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positioning  of 


the 


successful  procedure.  Moreover, 
measuring  volume  in  the  desired  location  is 


not  easy. 


Multi-Hole  Probe  System 


A  multi-hole  probe  is  rugged.  It  requires  only  one 
calibration  at  the  end  of  its  manufacturing  process.  Its 
use  is  very  simple  and  the  technique  can  provide 
information  on  all  the  primitive  variables  of  the  flow. 


Table  V.l  shows  the  comparison  matrix  of  each  main 
method  of  measurement  including  the  pitot-static  probe 
system . 


Table  V.l  Comparison  of  Different  Measurement  Technique 


Hot-Wire 

Anemometry 

Laser- 

Doppler 

velocimeter 

Multi  hole 
probe  system 

Pitot-static 
probe  + 
alpha  beta 
vane  system 

Calibration 

Often 

None 

Once 

Once 

Interference 

Low 

None 

Low/Medium 

Low/Medium 

Signal 

Continuous 

Discrete 

Continuous 

Continuous 

Frequency 

response 

High 

Medium 

Medium 

Medium/Low 

Difficulty 
of  operation 

Medium 

High 

Low 

Very  low 

Cost 

Medium 

High 

Low 

Very  low 

Accuracy 

High 

High 

High 

Medium/Low 

Use  on  FROG 

Almost 

impossible 

Impossible 

Very 

possible 

Very 

possible 
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1. 


5 -Hole  Probe 


A  5-hole  probe  was  custom  made  as  part  of  this  thesis 
search  for  better  air  data  sensors  onboard  the  FROG.  The  5- 
hole  probe  is  expected  to  be  more  precise  and  accurate  in 
measurement  of  a  -  (5  angles  compared  to  using  the  vanes 
mounted  on  the  pots.  Figure  V.7  shows  the  dimensional 
drawing  of  the  5-hole  probe. 
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Figure  V.7  Dimensional  drawing  of  5-hole  probe 

The  construction  of  probes  can  be  divided  into  the 

fabrication  of  internal  and  external  features.  The  external 
features  define  the  geometry  exposed  to  the  flow.  The 

internal  features  define  the  discreet  pressure  channels 
that  transmit  the  pressure  from  the  probe  tip  to  the 
pressure  transducers. 

a.  External  Features 

Figure  V.8  illustrates  the  arrangement  of  the 

probe  components.  Probe  tips  are  typically  made  of  brass  or 
stainless  steel.  There  are  three  parameters  defining  the 
tip  geometry  and  features:  tip  diameter,  shape  (conical, 
hemispherical,  and  faceted)  and  number  of  holes  (5  or  7)  . 
The  first  extension  (Figure  V.8)  precisely  matches  the 
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outside  diameter  of  the  tip  and  is  typically  up  to  20  tip 
diameters  long.  Ferrules  gradually  increase  the  outside 
diameter  of  the  probe  shaft,  providing  the  strength  of 
material  necessary  for  a  given  length.  If  needed,  a  second 
extension  can  be  added,  creating  a  longer  shaft.  A  third 
extension  can  also  be  added  if  needed,  creating  a  very  long 
probe  for  special  applications.  Mounts  are  typically 
hexagonal  to  allow  for  probe  roll  referencing. 


Figure  V.8  Typical  layout  of  5-hole  probe 


Figure  V.9  Details  of  5-hole  probe  tip 

b.  Internal  Features 

Each  of  the  holes  in  the  probe  tip  leads  to  a 
stainless  steel  tube,  with  its  inside  diameter  matching  the 
diameter  of  the  hole.  As  the  probe  shaft  diameter 
increases,  each  tube  is  telescoped  into  a  larger  tube. 
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which  finally  protrudes  from  the  back  of  the  mount.  Each 
connection  is  soldered  and  tested  for  strength  and  leakage. 
The  final  assembly  is  also  tested  for  pressure  "cross¬ 
talk",  i.e.  pneumatic  communication  between  two  or  more 
probe  holes  and  their  associated  tubing.  Figure  V.9 
presents  structural  details  of  a  typical  probe  tip  and 
sample  dimensions  for  a  0.062"  (1/16")  tip  diameter  probe 

tip . 


2.  Principle  of  Operation  for  Multi-Hole  Probe 

The  pressure  over  a  bluff  body  is  the  highest  at  the 
stagnation  point  and  lowest  near  separation.  If  the  flow 
direction  forms  small  angles  with  the  axis  of  the  probe 
(below  20°),  the  center  hole  registers  the  highest 

pressure.  If  however,  the  flow  is  steeply  inclined  with 
respect  to  the  probe,  as  shown  in  Figure  V.10  for  a  seven- 
hole  probe,  then  one  of  the  peripheral  holes  on  the 
windward  side  of  the  probe  tip  registers  the  highest 

pressure,  while  on  the  leeward  side  of  the  probe,  the  flow 

is  separated.  The  pressure  information  provided  by  the 
three  holes  in  the  separated  region  is  not  used.  Careful 

calibration  can  thus  allow  the  instrument  to  measure  the 
direction  and  the  magnitude  of  the  velocity. 
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Figure  V.10  Multi-hole  probe  inserted  in  flowfield 

Using  the  5-hole  probe  to  measure  a  and  /?  is  governed 
by  the  following  simplified  basic  principle:  Defining  the 

{P\  -  Pv  ) 

port  pressure  as  per  Figure  V.ll,  AP  is  defined  as  - ^ aig 

(. P2  +  P3  +  P4  +  P5 )  , 

where  Pa  = -  .  Therefore,  a  is  proportionate  to 


(P4-P5)  ,  _  .  .  (P2-P3)  ,  . 

-  whereas  p  is  proportionate  to  -  and  it  is 

A P  A P 

this  proportionality  that  allows  the  a  -  f5  to  be  determined 
once  the  calibration  is  carried  out. 


Figure  V.ll  Typical  5-hole  probe  port  numbering  system 
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Due  to  the  finite  size  of  the  probe  (and  the  very 
nature  by  which  it  operates) ,  one  may  anticipate  some 
interference  with  the  flow.  Caution  should  be  exercised  if 
one  attempts  to  measure  too  close  to  a  solid  wall  and  in 
boundary  layers.  Fluid  mechanics  phenomena  that  are  very 
sensitive  to  external  disturbances  have  also  been  known  to 
be  affected  by  the  probe.  For  example,  inserting  the  probe 
near  the  core  of  a  columnar  vortex,  like  the  tip  vortex  of 
a  wing,  may  induce  vortex  breakdown  and  inserting  the  probe 
in  a  laminar  shear  layer  or  a  laminar  boundary  layer  at 
moderate  Reynolds  numbers  may  induce  transition. 

3 .  Wind  Tunnel  Experiment 

A  wind  tunnel  experiment  was  conducted  to  compare  the 
a  -  (3  data  acquired  by  5-hole  probe-transducer  and  a  - /3 
vane-pots  systems.  Both  probes  were  mounted  in  the  wind 
tunnel  (which  has  a  flow  turbulence  level  of  0.23%)  to 
measure  only  /?  since  measuring  /?  alone  is  sufficient  for  a 
comparison  of  the  two  measurement  systems.  Figure  V.12 
shows  how  both  the  probes  were  mounted  in  the  wind  tunnel 
section.  The  free-stream  velocity  of  the  tunnel  was  set  to 
approximately  25  meters  per  second  since  the  FROG  operating 
speed  range  is  25  -  35  meters  per  second.  The  output 

signals  from  the  five  pressure  transducers  (Transducer  B 
type)  and  the  signals  from  the  /?  pot  were  processed  by  a 
processor  with  a  10-bit  analog  to  digital  converter  which 
was  sampling  data  at  a  freguency  of  40  hertz. 
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Figure  V.12  Set-up  of  /3  vane-pot  and  5-hole  probe- 

transducer  system 

Prior  to  the  experiment,  the  /?  pot  and  the  transducers 
were  calibrated  the  same  way  as  it  was  done  for  the  FROG  as 
mentioned  in  Chapter  II.  The  first  part  of  the  experiment 
involved  operating  the  wind  tunnel  at  a  specific  speed 
setting  and  measuring  /?  from  -8°  to  +12°  at  intervals  of  2° 
statically  while  the  second  part  involved  measuring  /?  from 
-10°  to  +12°  dynamically.  Since  five  transducer  B  types  were 
used  in  the  experiment,  the  port  pressures  (PI,  P2,  P3,  P4 

and  P5)  were  all  referenced  to  atmospheric  pressure.  The 
pressure  values  were  reduced  to  compute  the  /?  values  using 
a  calibration  software.  Multiprobe®.  Note  that  Multiprobe® 
produces  many  parameters  of  interest:  a  -  /3 ,  u-v-w 

components  of  velocity,  total  velocity  V,  total  pressure 
and  static  pressure  and  total  temperature  and  static 
temperature . 
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The  (3  vane-pot  system  was  mounted  on  a  turn-table  as 
shown  in  Figure  V.12.  The  zero  /?  reference  was  obtained  by 
aligning  the  probe  to  be  parallel  to  the  wind  tunnel  and 
the  P  vane  to  be  parallel  to  the  airstream.  The  probe 
alignment  was  done  visually;  standing  in  the  middle  of  the 
wind  tunnel  and  adjust  the  turn-table  till  the  probe  is 
aligned.  Once  the  probe  was  aligned,  the  vane  was  aligned 
by  running  the  wind  tunnel  up  to  speed  and  let  the 
airstream  align  the  vane.  Once  that  was  done,  the  turn¬ 
table  angular  scale  was  reset  to  zero  as  reference. 
Subseguently,  as  P  was  measured  in  intervals  of  2°,  it  was 
based  on  the  turn-table  angular  scale. 

4 .  Wind  Tunnel  Experiment  Results 

Results  obtained  from  the  first  part  of  the  experiment 
are  as  shown  in  Figure  V.13.  From  Figure  V.13,  it  can  be 
assessed  that  both  systems  of  P  measurement  were  in 
reasonably  close  agreement  although  there  were  slight 
differences  between  the  P  values  especially  at  higher  P 
value . 
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Figure  V.13  Static  measurement  of  /?  from  5-hole  probe  and 

P  pot 

Table  V.2  summarizes  the  mean  and  RMS  values  of  the 
signal  from  the  /?  vane-pot  and  5-hole  probe-transducer 
system.  Since  the  ft  vane-pot  readings  were  used  as 
reference  (the  visual  alignment  in  the  wind  tunnel  was 
taken  as  zero  for  the  beta  vane-pot  system)  ,  the  mean 
values  were  corrected  so  that  a  better  comparison  can  be 
made.  Table  V.3  shows  the  corrected  mean  and  difference 
values  as  well  as  the  RMS  values  of  the  signals  from  both 
systems . 
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Table  V.2  Mean  and  RMS  values  from  /?  vane  and  5-hole  probe 


Mean  value  from 
(3  vane 

Mean  value  from 
5-hole  probe 

RMS  value  from  (5 

vane 

RMS  value  from  5- 
hole  probe 

o 
\ — 1 

CO 

1 

-8 . 12 

0.29 

0.38 

CM 

O 

CD 

1 

-5.82 

0.29 

0.38 

1 

o 

CD 

-3.89 

0.26 

0.42 

-1 . 94 

-1.56 

0.34 

0.39 

0.06 

0.42 

0.29 

0.31 

1 . 86 

2.39 

0.29 

0.33 

3.91 

4.31 

0.28 

0.29 

5.96 

6.46 

0.29 

0.36 

8.01 

8.67 

0.27 

0.37 

9.94 

10.69 

0.30 

0.36 

12.01 


12 . 61 


0.27 


0.35 


Table  V.3  Corrected  Mean  and  RMS  values  from  /?  vane  and  5- 
hole  probe 


Tunnel 

reference 

Corrected  mean 
value  from  /? 
vane 

Corrected  mean 
value  from  5- 
hole  probe 

Difference 

RMS  value 
from  P 
vane 

RMS  value 
from  5-hole 
probe 

O 

O 

CO 

1 

-8.16 

-8.18 

-0.02 

0.29 

0.38 

o 

o 

to 

1 

-6.08 

-5.88 

0.20 

0.29 

0.38 

-4.00 

-4 . 15 

-3.95 

0.21 

0.26 

0.42 

o 

o 

CM 

1 

-2.00 

-1 . 62 

0.38 

0.34 

0.39 

0.00 

0.00 

0.36 

0.36 

0.29 

0.30 

2.00 

1 . 80 

2.33 

0.53 

0.29 

0.33 

4.00 

3.85 

4.25 

0.40 

0.28 

0.29 

6.00 

5.90 

6.40 

0.50 

0.29 

0.36 

8.00 

7 . 95 

8.61 

0.67 

0.27 

0.37 

10.00 

9.89 

10.64 

0.75 

0.30 

0.36 

12.00 

11 . 95 

12.55 

0.60 

0.27 

0.35 

Note  that  the  RMS  values  of  the  5-hole  probe- 

transducer  system  are  larger  than  those  of  the  P  vane-pot 
system.  This  suggests  that  in  terms  of  precision,  the  /? 
vane-pot  system  has  better  performance.  In  an  effort  to 
explore  the  possibility  of  a  scaling  relationship  between 
the  P  values  from  the  P  vane-pot  and  5-hole  probe- 

transducer  system,  the  mean  values  were  plotted  against 
each  other  as  shown  in  Figure  V.14.  A  curve  fitted  eguation 
of  P  (from  P  vane)  =  0.9692  x  /?  (from  5-hole  probe) 

0.3456  was  obtained  and  is  not  unreasonably  far  from  the 
ideal  equation  of  P  (from  /?  vane)  =  1.0  x  /?  (from  5-hole 

probe)  . 
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Figure  V.14  Relationship  of  /?  values  from  (3  vane-pot  and 
5-hole  probe-transducer  system 

The  scaling  relationship  was  applied  to  the  /?  signal 
in  Figure  V.12.  Figure  V.14  shows  the  results.  Note  that 
both  signals  were  in  very  close  agreement  after  the  scaling 
was  applied  to  the  signal  from  the  5-hole  probe-transducer 
system . 


Figure  V.15  Scaled  data  from  5-hole  probe 
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A  PSD  analysis  was  conducted  and  the  PSD  plots  are  as 
shown  in  Figure  V.16.  For  comparison  purposes,  the  PSD  of 
the  test-flight  [3  signal  from  the  (3  vane-pot  system  on  the 
FROG  was  superimposed  (green  curve  in  Figure  V.16) .  Note 
the  significant  larger  power  of  the  test-flight  (3  signal  as 
compared  to  those  from  the  [3  vane-pot  and  5-hole  probe- 
transducers  in  the  wind  tunnel  experiment.  This  can  be 
accounted  by  the  fact  that  aircraft  vibrations  and  other 
flight  operating  environment  effects  were  present  during 
the  test  flight  (which  is  not  present  in  the  wind  tunnel 
experiment)  .  Figure  V.16  also  shows  that  both  the  5-hole 
probe-transducer  system  and  /?  vane-pot  system  have 
approximately  the  same  bandwidth.  The  /?  vane-pot  system  has 
better  roll-off  performance  than  the  5-hole  probe- 
transducer  system.  In  terms  of  out-of-band  frequencies 
(approximately  above  2  hertz) ,  the  5-hole  probe-transducer 
has  better  attenuation  than  the  [3  vane-pot  system. 


Figure  V.16  PSD  plots  of  f3  signal  from  5-hole  probe  and  f3 
pot  (wind  tunnel  and  on  the  FROG) 
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Figure  V.17  shows  the  result  of  the  second  part  of  the 
experiment  which  was  dynamically  measuring  /?  by  sweeping 
from  -10°  to  12°  and  back  to  -10°.  Note  that  the  sweep  of  (3 
from  -10°  to  12°  occurred  at  a  different  sweep  rate  from  the 
sweep  of  [3  from  12°  back  to  -10°  as  the  turntable  mount  of 
the  wind  tunnel  has  a  sweep  rate  of  4.5°  per  second  when 
1 3  is  increasing  and  5°  per  second  when  [3  is  decreasing. 


Figure  V.17  Dynamic  measurement  of  [3  from  5-hole  probe  and 

(3  pot 

Upon  careful  scrutiny  of  Figure  V.17,  as  [3  was  swept 
increasing  from  -10°  to  12°  at  a  rate  of  4.5°  per  second, 
the  signal  from  the  5-hole  probe-transducer  system  was 
leading  the  signal  from  the  (3  vane-pot  system  for  most  of 
the  time;  this  can  be  observed  from  the  400th  to  the  950th 
data  point.  When  (3  was  swept  decreasing  from  12°  back  to  - 
10°  at  a  rate  of  5°  per  second,  the  signal  from  the  [3  vane- 
pot  system  was  instead  leading  the  signal  from  the  5-hole 
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probe-transducer  system;  note  the  950th  to  the  1200th  data 
point.  A  hysteresis  effect  is  observed  here.  In  this 
dynamic  case,  a  parameter  of  interest  for  the  5-hole  probe- 
transducer  system  will  be  the  lag  constant  which,  in  [1]  is 
given  as 


X  = 


12%/uLV 

nd4P 


(12) 


Where : 


//  is  the  coefficient  of  viscosity  of  air 
L  is  the  total  length  of  tubing  used 

V  is  the  combined  volume  in  tubing  and  transducer 
cavity 

D  is  the  internal  diameter  of  the  tubing  used 

P  is  the  reference  pressure  in  the  tubing 

Since  a  total  length  of  approximately  6  feet  of  tubing 
(internal  diameter  1/32  inch  was  used),  the  lag  constant 
for  the  5-hole  probe-transducer  was  determined  to  be  0.12 
seconds.  For  applications  of  pressure  measurements,  it  is 
reasonable  to  assume  that  the  lag  constant  of  0.12  seconds 
is  not  an  issue. 

5.  Use  of  5-Hole  Probe  on  the  FROG 

The  pressure  configuration  on  the  5-hole  probe  used  in 

the  wind  tunnel  experiment  was  that  each  of  the  port 
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pressures  was  referenced  to  the  atmospheric  pressure.  For 
the  FROG,  a  few  configurations  are  possible.  The  first 
configuration  is  such  that  five  absolute  pressure 
transducers  are  used  for  measuring  the  port  pressures. 
Another  possible  configuration  will  be  using  transducer  B 
type  (differential  type)  and  referencing  each  of  the  port 
pressures  to  the  static  pressure  of  the  FROG;  however,  an 
absolute  pressure  value  of  the  static  pressure  is  required 
(which  implies  the  use  of  five  differential  transducers  and 
one  absolute  transducer)  by  Multiprobe®  for  data  reduction. 

The  last  possible  configuration  is  using  only  four 
transducer  B  type  and  referencing  all  the  port  pressures  to 
the  center  port  pressure.  This  configuration  allows 
Multiprobe®  to  compute  a  -  (5  from  only  the  u-v-w  component 
of  velocity.  But  the  total  pressure  and  static  pressure 
values  reduced  from  Multiprobe®  cannot  be  used  as  the 
reference  pressure  in  this  case  will  be  unknown.  But  since 
the  subject  of  interest  is  really  only  a  -  (3 ,  this 
configuration  is  recommended  for  use  on  the  FROG  as  it 
minimized  the  number  of  transducers  used,  allows  simpler 
pressure  tubings  installation  and  reduces  the  number  of 
channels  for  data  acquisition. 
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VI  DISCUSSION 


A.  SENSOR  CALIBRATION 

The  calibration  of  the  pressure  transducer  and  the  a- 
(3  pots  is  of  great  importance,  as  a  good  calibration  will 
give  accurate  air  data  that  is  desired.  A  good  engineering 
practice  is  that  the  sensors  are  to  be  calibrated  in  the 
laboratory  (just  the  sensors  with  the  same  processor 
acquiring  data  in  flight)  and  on  the  aircraft  before  and 
after  test  flights.  The  three  calibration  results  are  then 
evaluated,  compared  and  modified  as  necessary.  As  far  as 
this  thesis  is  concerned,  the  sensors  were  calibrated  in 
the  laboratory  and  on  the  FROG  before  the  test  flight.  The 
calibration  results  were  in  close  agreement.  However,  as 
the  FROG  was  lost  due  to  an  air  accident  during  the  test 
flight  on  9  Oct  2001,  the  calibration  after  the  test  flight 
could  not  be  carried  out.  Nevertheless,  the  initial  two 
calibrations  carried  out  were  deemed  to  be  sufficient. 

As  far  as  the  calibration  of  the  5-hole  probe  is 
concerned,  it  is  not  an  issue  as  it  is  done  by  the 
manufacturer.  The  calibration  for  the  5-hole  probe  is  a 
highly-accurate,  more  than  two  thousand  point  calibration. 
It  also  includes  a  post  calibration  error  analysis  and 
quality  control.  Since  the  FROG  is  operating  at  the  low 
Mach  number  range,  calibration  at  a  single  speed  is 
sufficient  for  the  5-hole  probe. 
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B.  RESOLUTION 

Since  the  air  data  signal  had  been  digitized,  one 

issue  that  arises  is  the  resolution  or  quantization  level 
of  the  digitized  signal.  The  digitized  signal  is  an  integer 
number  of  counts  basically.  The  resolution  of  a  signal  is 
exactly  one  count  which  corresponds  to  an  engineering  unit 
that  depends  on  the  sensor  calibration.  In  the  case  of  the 
FROG,  the  engineering  units  of  interest  are  meters  per 

seconds  (airspeed)  and  degrees  (a-/3)  .  As  mentioned  in 
Chapter  II,  the  resolution  of  the  pressure  transducer  is 
approximately  0.008  meters  per  second  and  0.09°  for  the  a  - 
/?  pots.  It  was  concluded  from  [6]  that  as  long  as  the 

resolution  magnitude  is  much  smaller  than  the  noise  level 
of  the  digitized  data,  it  is  very  much  acceptable  since  it 
allows  detection  of  small  variations  of  the  analog  input 
signal . 

The  data  used  for  plotting  Figures  III. 5  through 
III. 10  which  presented  the  air  data  signal  before  and  after 
filtering,  was  used  to  determine  the  average  nose  level. 
The  average  noise  level  for  the  pressure  transducer 

(transducer  A)  was  found  to  be  26  digital  counts  which 
translate  to  a  1.25  meters  per  second  noise  level.  The 
average  noise  level  of  for  the  a  -  (3  pots  was  found  to  be 
0.1°.  It  is  obvious  that  the  resolution  magnitude  of  the 
pressure  transducer  is  significantly  lesser  than  its  noise 
level.  Therefore,  as  far  as  the  resolution  level  for  the 
pressure  transducer,  it  is  much  acceptable  based  on  the 
processor  with  a  12-bit  analog  to  digital  converter.  But 
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the  a  -  [3  pots  do  require  an  analog  to  digital  converter 
that  is  more  than  12-bit. 

C.  REAL-TIME  PROCESSING  OF  AIR  DATA 

The  air  data  in  this  project  was  post-processed  rather 
than  real-time  processed.  The  fact  that  air  data  was  post- 
processed  presented  a  great  degree  of  simplicity  for  this 
project  as  far  as  processing  of  air  data  is  concerned.  If 
real-time  processing  of  air  data  is  desired,  MATLAB®  Real 
Time  Workshop  can  be  used.  The  research  effort  in  real-time 
processing  of  air  data  is  currently  being  pursued  in  [9]  . 
In  Real  Time  Workshop,  filtering,  re-sampling,  calibration 
and  corrections  sub-modules  have  to  be  created  in  order  for 
it  to  perform  real-time  processing.  In  particular,  the 
calibration,  static  and  dynamic  correction  results  from 
this  thesis  have  to  be  implemented  and  created  as  a  sub- 
module  in  Real  Time  Workshop.  Filtering  and  re-sampling 
also  have  to  be  built  in  as  a  sub-module  into  Real  Time 
Workshop.  Real-time  processing  of  air  data  also  can  be 
integrated  together  with  the  wind  estimation  model  as 
mentioned  in  Chapter  III  Section  F  to  provide  real  time 
wind  estimation  which  is  crucial  for  Chem/Bio  attack 
response.  Real-time  processing  of  air  data  with  proper 
display  provides  the  pilot  of  the  FROG  with  a  form  of  real 
time  feedback  which  is  also  beneficial. 

D.  PRESSURE  TRANSDUCER  A  VERSUS  PRESSURE  TRANSDUCER  B 

Chapter  V  covered  in  detail  the  experiment  conducted 
to  compare  both  transducers  A  and  B.  Both  transducers 
exhibit  the  same  bandwidth  and  roll-off.  But  in  terms  of 
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out-of-band  frequencies,  transducer  B  does  suggest  that  is 
has  better  attenuation.  However,  to  be  complete,  issues 
pertaining  to  other  transducers  characteristics  like  age, 
linearity,  repeatability  and  environmental  effects  have  to 
be  also  thoroughly  researched.  In  terms  of  frequency 
response  and  noise  characteristics,  the  wind  tunnel 
experiment  results  as  presented  in  Chapter  V  do  suggest 
that  transducer  B  is  a  more  favorable  choice  for  future  use 
in  the  FROG. 

E.  a-/3  VANE-POT  SYSTEM  VERSUS  5-HOLE  PROBE -TRANSDUCER 

SYSTEM 

The  details  of  the  results  from  the  wind  tunnel 
experiment  conducted  to  compare  both  a  -  (3  vane-pot  and  5- 
hole  probe-transducer  systems  were  covered  in  Chapter  V.  As 
mentioned  before,  the  RMS  value  of  the  signal  achieved  by 
the  5-hole  probe-transducer  system  is  larger  than  that  of 
the  P  vane-pot  system.  Since  a  lower  RMS  value  is  desired, 
this  indicated  that  the  /?  vane-pot  system  is  better  in  this 
aspect.  For  frequency  response,  both  the  5-hole  probe- 
transducer  system  and  /?  vane-pot  system  has  approximately 
the  same  bandwidth.  The  /?  vane-pot  system  has  better  roll¬ 
off  performance  than  the  5-hole  probe-transducer  system  and 
in  terms  of  out-of-band  frequencies;  the  5-hole  probe- 
transducer  has  better  attenuation  than  the  /?  vane-pot 
system . 

Both  the  measurement  systems  have  to  be  subjected  to  a 
real  flight  operating  environment  and  be  tested  in  flight. 
This  will  provide  a  more  comprehensive  review  and 
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comparison  of  both  measurement  systems  as  the  experiment  in 
the  wind  tunnel  alone  is  not  sufficient  to  simulate  the 
real  flight  operating  environment,  for  example,  the  noise 
level  in  wind  tunnel  and  in  flight  definitely  differs  a 
significant  amount  for  the  a  -  (3  vane-pots  system.  However, 
the  wind  tunnel  experiment  did  reveal  an  important  fact: 
the  5-hole  probe-transducer  is  capable  of  measuring  a  -  f3  on 
the  FROG. 
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VII  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSION 

The  air  data  system  using  standard  sensors  was 
installed  on  the  FROG  after  consideration  of  air  data  error 
associated  with  several  locations  of  the  sensors  on  the 
FROG.  A  test  flight  was  carried  out  and  air  data  was 
acquired  successfully.  The  post  processing  of  the  air  data 
which  include  filtering,  re-sampling,  calibrating,  static 
and  dynamic  correction  (numerical  simulation  for  static 
source  correction)  were  carried  out  as  well.  Preliminary 
wind  extraction  was  carried  out  (refer  to  [9])  from  the 
post  processed  air  data  and  results  are  very  promising  as 
the  preliminary  wind  estimation  was  found  to  be  in  close 
agreement  with  wind  data  from  the  Meteorology  Department. 
Further  wind  extraction  will  be  carried  out  and  researched 
in  [9]  . 

Pressure  transducers  were  compared  by  running  wind 
tunnel  experiments  and  analyzing  the  quality  of  output 
signals  from  the  transducers. 

A  5-hole  probe-transducer  system  was  also  used  in  a 
wind  tunnel  experiment  to  compare  with  the  a  -  ft  vane-pot 
system.  The  performance  of  the  5-hole  probe-transducer  and 
the  a  -  ft  vane-pot  system  were  analyzed  and  discussed  in 
terms  of  RMS  values  and  frequency  response. 

Specific  issues  pertaining  to  air  data  acquisition 
were  also  explored  and  discussed.  This  thesis  project  will 

pave  the  way  for  the  UAV  to  be  used  for  Chem/Bio  attack 
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response.  Use  of  a  UAV  to  make  in  situ  measurements  of 
particulates  as  well  as  air  and  inertial  data  for  winds 
extraction  is  a  unique  contribution  to  the  area  of  homeland 
security . 

B.  RECOMMENDATIONS 

With  the  success  completion  of  this  thesis  effort,  the 
following  can  be  pursued  in  the  near  future: 

1.  The  change-out  and  addition  of  specific 
components  within  the  air  data  system;  a  16-bit  or  higher 
analog  to  digital  converter  in  the  processor  can  be  used 
instead  of  the  current  12-bit  analog  to  digital  converter. 
Since  pressure  transducer  B  has  better  frequency  response 
and  is  less  noisy  than  transducer  A,  it  is  also  recommended 
that  pressure  transducer  B  be  used  for  the  pitot  static 
pressure  measurement  and/or  5-hole  probe  pressure 
measurement.  Since  the  density  correction  of  the  air  data 
was  carried  out  with  data  from  the  Meteorology  Department 
(as  mentioned  in  Chapter  III),  it  is  also  recommended  that 
the  FROG  be  installed  with  pressure  altitude  and 
temperature  sensors  so  that  the  exact  density  of  the  FROG 
at  a  particular  altitude  can  be  computed.  This  allows  a 
more  accurate  density  correction  to  be  carried  out. 

2 .  Use  of  more  advanced  numerical  software  for  the 
dynamic  correction;  PSW®  was  used  to  make  off-body  flow 
analysis  to  calculate  the  corrections  to  be  made  to  the  air 
data.  The  software  suite  had  since  been  updated  to  a  newer 
windows  version.  Therefore,  it  is  recommended  that  the 
newer  version  or  more  advanced  numerical  software  be  used. 
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The  challenge  for  more  advanced  numerical  software  however 
will  be  to  re-create  the  numerical  model  (grid)  of  the  FROG 
in  the  new  software. 

3.  Using  Real  Time  Workshop  in  MATLAB®,  establish  a 
real-time  air  data  processing  system  instead  of  the  current 
way  of  post  processing  air  data.  It  is  also  recommended 
that  Multiprobe®  be  integrated  into  Real  Time  Workshop  in 
MATLAB®  to  allow  real-time  data  reduction  from  the  pressure 
readings  obtained  from  the  5-hole  probe.  This  requires 
consultation  with  the  manufacturer  of  the  5-hole  probe. 

4.  Test  flight  both  the  /?  vane-pot  and  5-hole  probe- 
transducer  on  the  FROG  to  compare  the  performance  of  both 
a  - measurement  systems.  Together  with  the  findings  from 
the  wind  tunnel  experiment  conducted  in  this  project,  a 
decision  can  be  then  made  to  select  the  appropriate  a- ft 
measurement  system. 
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APPENDIX  A.  DESCRIPTION  OF  FOG-R  UAV  EXTRACTED  FROM 

[10] 


Tactically  Expendable 
Remote  Navigator 
(TERN) 

TERN  is  an  acronym  for  "Tactically 
Expendable  Remote  Navigator.  This 
versatile  UAV  is  capable  of  a  variety 
of  missions,  recoverable  in  peacetime  training  missions,  yet  of  sufficiently  low 
cost  to  be  discarded  after  actual  battle  if  necessary.  A  variant  of  the  vehicle, 
FOG-R,  has  been  flown  non-line-of-sight  (NLOS)  "over  and  beyond  the  hill,"  using 
a  Fiber  Optic  Data  Link  (FODL)  for  command  uplink  and  video  downlink.  This 
FOG-R  version  is  jam-proof.  Another  variant,  TERN-C,  has  carried  a  remote 
sensor  for  detecting  chemical  warfare  agents. 


The  TERN  is  designed  to  carry  up  to  22  lbs. 

(10  kg)  of  payload  for  periods  of  up  to  four 
hours.  The  wing  is  fitted  with  flaps  which  can 
be  trimmed  to  provide  slow  flight  speed  for 
surveillance  and  extended  to  facilitate  tight 
landings. 

TERN's  optional  equipment  includes  a 
programmable  autopilot,  GPS  navigation,  video 
system,  and  other  electronic  sensors. 
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APPENDIX  B.  DETAILED  SPECIFICATIONS  OF  SENSOR 
EXTRACTED  FROM  [11] 


VISHAY 

T 


1/2”  (12.7mm)  Single  ■  Turn  Wirewound 
Precision  Potentiometer 


Model  140, 142 

Vishay  Spectrol 


Model  140  Model  142 

Bushing  Mount  Servo  Mount 


ELECTRICAL  SPECIFICATIONS 


FEATURES 

•  50ftto20Kfl 

•  Bushing  or  Servo 


Total  Resistance 

STANDARD 

50ft  to  20  KQ 

Tolerance 

±5% 

Absolute  Minimum  Resistance 

Linearity  x  Total  Resistance  or 

0.5ft,  whichever  is  greater 

Linearity  (Independent) 

±1.0% 

Noise 


Power  Rating 


2  watts  at  40°C  ambient  derating  linearly  to  zero  at  125°C 


Insulation  Resistance 

I.OOOMft  min.  0.500VDC 

Dielectric  Strength 

1.000V rms.  60  Hz 

Rotation 

ua 

320°  ±  5® 

142 

350°  +  0°  -  4° 

End  Voltage 

Linearity  x  total  applied  voltage  for  total  resistance  above  20ft; 

2.0%  of  total  applied  voltage  for  20ft  and  below 

MATERIAL  SPECIFICATIONS 


Shaft 

Stainless  steel,  non  magnetic 
non-passivated 

Housing 

Aluminum,  anodized 

Rear  Lid 

Molded  glass  filled 
thermoset  plastic 

Terminals 

Brass,  gold  plated 

Mounting  Hardware 
Lockwasher  Internal  Tooth: 
Panel  nut: 

(Model  140  only) 

Steel,  nickel  plated. 

Brass,  nickel  plated 

ENVIRONMENTAL  SPECIFICATIONS 


20G  thru  2000  Hz 


Vibration 


Shock 


Salt  Spray 
Rotational  Life 


Temperature  Range 


50g 


%  Hours 


500,000 

Shaft  Revolutions 


ORDERING  INFORMATION _ _ 

This  part  number  consists  of  four  groups  of  digits.  The  first  group  is  the  Spectrol  model  number.  The  second  digit  describes  the  mechani¬ 
cal  options  available.  The  third  digit  describes  the  other  optional  features  which  can  be  supplied.  The  fourth  group  is  the  standard  EIA 
resistance  code. 

Example:  140  -  0  -  0  -  203 

140  0  0  203 

MODEL  MECHANICAL  OPTIONS  OTHER  OPTIONAL  RESISTANCE  CODE 

FEATURES 


0.  Stops,  Slotted  Shaft  (std) 

1.  Plain  Shaft 

2.  Shaft  Lock 

3.  Continuous  Rotation 

4.  Combination  1  &  2 

5.  Combination  1  &  3 

6.  Combination  2  &  3 

7.  Combination  1 , 2,  &  3 

0.  Continuous  Rotation. 
Plain  Shaft  (std.) 


0.  Standard  Torque 

1 .  Center  Tap  (10K  max.  Rt) 

2.  High  Torque 

3.  Sealed  Construction 

4.  Combination  1  &  2 

5.  Combination  1  &  3 

6.  Combination  2  &  3 

7.  Combination  1 . 2,  &  3 

0.  Standard  Torque 
1 .  Center  Tap  (10K  max.  Rt) 


2  -  First  Significant  figure 

0  =  Second  significant  figure 

3  =  Number  of  zeros  following 


Same  as  above 


Example:  Part  number  140-0-0-203  describes  bushing  mount  140  with  stops,  slotted  shaft,  standard  torque  and  resistance  of  20KQ 
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Model  140,142 _ 

Vishay  Spectral  1/2”  (1 2.7mm)  Single  -  Turn  Wirewound 

Precision  Potentiometer 


VISHAY 


DIMENSIONS  in  inches  (millimeters) 


*0  0000 


MECHANICAL  SPECIFICATIONS 

PARAMETER 

Rotation 

140 

142 

330°  ±  5° 

360°  continuous 

Bearing  Type 

SLEEVE  BEARING 

BALL  BEARING 

Torque  (Maximums) 

Starling 

0.2  02  -  in  (14.40gm  -  cm) 

0.075  oz  -  in  (5.40gm  -  cm) 

Running 

0.2  oz  •  in  (14.40gm  -  cm) 

0.05  oz  -  in  (3.60gm  •  cm) 

Traversing 

- 

- 

Dead  Zone 

Not  applicable 

0.20  oz  -  in  (14.40gm  •  cm) 

Stop  Strength 

5  in  -  lbs  (5.76  kgm  •  cm)  static  (140  only) 

Runouts  (Maximum) 

140 

142 

Shatt  (TIR) 

0.002  in  (0.05cm) 

0.002  in  (0.05cm) 

Pilot  Dia  (TIR) 

0.002  in  (0.05cm) 

0.002  in  (0.05cm) 

Lateral  (TIR) 

0.003  in  (0.08cm) 

0.002  in  (0.05cm) 

Shatt  End  Play 

0.006  in  (0.15cm) 

0.004  in  (0.10cm) 

Shaft  Radial  Play 

0.003  in  (0.08cm) 

0.002  in  (0.05cm) 

MAXIMUM 

MAXIMUM 

STD 

CURRENT 

VOLTAGE 

RESISTANCE 

RESO- 

OHMS 

AT  40°C 

ACROSS 

VALUES 

LUTION 

PER 

AMBIENT 

COIL 

COEF. 

w 

(%) 

TURN 

tmA) 

(V) 

(ppm/'C) 

50 

0.542 

0.271 

200.0 

10.0 

20 

100 

0.431 

0.431 

141.0 

14.1 

20 

0.361 

0.722 

100.0 

20.0 

20 

0.312 

1.56 

63.2 

31.6 

20 

0.255 

2.55 

44.7 

44.7 

20 

0.197 

3.94 

31.6 

632 

20 

5K 

0.170 

8.50 

20.0 

100.0 

20 

0.147 

14.7 

14.1 

141.0 

20 

0.105 

21.0 

10.0 

200.0 

20 
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APPENDIX  C.  OVERVIEW  OF  PERSONAL  SIMULATION  WORKS 
EXTRACTED  FROM  [4]  (LOFTSMAN,  CMARC  AND  POSTMARC) 


A.  GENERAL 

Personal  Simulation  Works  is  a  PC  based  software  suite 
that  provides  for  the  three  primary  CFD  requirements;  3D 
modeling  of  an  aircraft  (LOFTSMAN) ,  panel  code  flow  solver 
(CMARC)  ,  and  post-processing  of  the  computed  flow  field 
(POSTMARC) .  The  software  package  contains  three 
applications  hosted  on  the  IBM  compatible  personal 
computer.  Each  software  program  is  discussed  separately. 

B .  LOFTSMAN 

LOFTSMAN  is  a  Windows  hosted  ad  modeling  tool  that 
generates  surface  panel  distributions  for  CMARC  or  PMARC 
input  files.  The  program  is  based  on  conics,  which  allows 
rapid  lofting  of  streamlined  bodies  such  as  aircraft 
fuselages  and  engine  nacelles.  In  addition,  wing  and 
control  surfaces  can  be  desired  with  the  extensive  library 
of  airfoil  templates  or  with  user  specified  coordinates. 
The  software  is  well  decremented,  including  a  tutorial,  in 
the  Personal  Simulation  Works  User  Guide.  LOFTSMAN  is 
primarily  designed  for  creating  new  objects,  but  an 
existing  airframe  can  be  matched  quite  closely  with  just  a 
detailed  three-view  drawing  that  includes  frame  cross 
sections . 
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1. 


Streamlined  Bodies 


LOFTSMAN  functionality  is  divided  into  Body  Objects 
and  Wing  Objects.  In  general,  they  remain  separate  unless 
the  intersection  between  a  wing  and  body  is  required.  Body 
Objects  are  created  using  a  family  of  curves  called  second- 
degree  conics.  Circles,  ellipses,  parabolas  and  hyperbolas 
are  among  this  group.  An  entire  fuselage  is  described  by 
specifying  just  four  lines.  These  are  the  top  waterline 
(TW)  ,  bottom  waterline  (BW) ,  the  maximum  breadth  line  (MB) 
and  the  waterline  of  the  maximum  breadth  line  (WW)  .  For 
each  line,  the  beginning,  ending  and  a  few  points  along  the 
line  are  specified.  Control  points  are  also  specified  with 
a  curvature  factor  that  allows  LOFTSMAN  to  generate  a 
smooth  conic  between  the  points. 

2 .  Wings  and  Control  Surfaces 

Wings  and  control  surfaces  are  easily  specified  in 
LOFTSMAN  using  a  short  input  file  created  with  any  text 
editor.  The  file  specifies  root,  intermediate  and  tip  rib 
section,  location,  axis,  chord  and  incidence.  LOFTSMAN  then 
fairs  a  smooth  surface  through  the  rib  sections.  Washout  is 
specified  by  varying  the  incidence  of  the  root  and  tip 
ribs.  Sweep-back  is  controlled  by  staggering  the  tip  rib 
location  with  respect  to  the  root  ribs.  Once  the  general 
wing  surface  is  specified,  control  surfaces  such  as 
ailerons,  flaps  and  elevators  can  be  deflected  and  meshed. 
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3. 


Patches 


LOFTSMAN  automatically  meshes  3D  surfaces  and  creates 
patches  for  CMARC/PMARC  input  files.  The  distinction 
between  a  mesh  and  a  patch  is  important.  A  mesh  is  a  set  of 
quadrilateral  and  triangular  panels  that  represent  the 
surface  of  a  wing  or  body.  When  the  set  of  panels  is 
organized  and  formatted  to  create  a  sub-component  portion 
of  a  CMARC  or  PMARC  input  file,  it  is  called  a  patch. 

A  body  or  wing  surface  is  first  meshed  at  a  density 
specified  by  the  user.  Panel  compression  options  include 
cosine  and  half-cosine  spacing.  After  meshing  the  object, 
one  saves  it  to  a  text  file  as  a  formatted  patch.  One  then 
opens  the  patch  file  with  any  text  editor  and  copies/pastes 
the  patch  text  into  the  appropriate  location  in  the  CMARC 
input  file. 

Each  control  surface  deflection  requires  a  separate 
mesh  and  formatted  patch.  For  instance,  to  evaluate  roll 
performance  one  needs  to  separately  mesh  an  upward  aileron 
deflection  on  the  right  wing  and  a  downward  deflection  on 
the  left  wing.  If  multiple  deflections  of  a  single  control 
surface  are  required,  each  deflection  must  be  meshed 
separately . 

C.  CMARC 

CMARC  is  the  C  version  of  PMARC  low-order  panel  code. 
Inviscid,  irrotational ,  incompressible,  potential  flow  is 
assumed.  Low-order  means  that  source  and  doublet  strength 
distribution  is  constant  across  a  panel.  There  is  no 
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attempt  to  match  that  source  or  doublet  strength  of  an 
adjacent  panel  at  a  common  edge.  Advanced  features  include 
internal  flow  modeling  and  time  stepping  wake  models. 

PMARC  version  12.19  was  released  as  FORTRAN  77  source 
code  in  1992.  CMARC  was  rewritten  in  the  C  language  and 
compiled  for  hosting  on  IBM  compatible  personal  computers 
by  Aerologic,  Inc.  The  program  runs  under  the  DOS  operating 
system.  It  will  also  run  in  a  DOS  window  from  Windows  3.1, 
95,  98,  2000  or  XP .  Enhanced  features  include  command-line 
options  and  flexible  memory  management.  Command  line 
options  simplify  batch  processing  by  adding  an  extensive 
set  of  switches  that  can  be  set  external  to  the  CMARC  input 
ale.  Flexible  memory  management  provides  for  the  automatic 
sizing  of  arrays  without  having  to  recompile  the  source 
code . 

D .  POSTMARC 

POSTMARC  is  a  Windows  post-processing  program  for  the 
visualization  of  CMARC  and  PMARC  output  files.  Capabilities 
include  body  geometry,  wake  stepping,  surface  pressure  and 
streamline  visualization.  POSTMARC  also  provides  the 
capability  to  integrate  pressure  and  skin  action  forces 
over  the  model  geometry.  This  proves  particularly  useful 
when  one  desires  to  recalculate  loads  around  a  different 
center  of  gravity. 

An  interesting  feature  for  design  work  is  the 
integration  of  panel  surface  area  to  total  wetted  area. 
After  lofting  a  new  geometry  in  LOFTSMAN,  a  quick  check  of 
geometry  is  made  by  running  CMARC  with  the  in  command  line 
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toggle.  The  total  wetted  area  is  then  checked  in  POSTMARC. 
This  function  is  particularly  useful  when  working  to  reduce 
skin  friction  drag. 

Versions  1.17.3  and  later  of  POSTMARC  include  the 
capability  to  integrate  skin  friction  drag  coefficient  over 
the  model  geometry.  It  is  important  to  note  that  a  key 
piece  of  the  drag  equation  is  missing  from  a  POSTMARC 
solution.  CMARC  provides  induced  drag  from  the  surface 
pressure  distribution  and  skin  friction  drag  from  the  2D 
boundary  layer  code.  Skin  friction  is  only  calculated  up  to 
the  point  of  boundary  layer  separation.  Pressure  drag  due 
to  separation,  a  major  portion  of  the  drag  equation,  is 
missing  from  a  CMARC /POSTMARC  solution. 
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APPENDIX  D.  DETAILED  SPECIFICATIONS  OF  PRESSURE 
TRANSDUCER  B  EXTRACTED  FROM  [12] 


SERIES  850  BOARD-MOUNTABLE  PRESSURE  TRANSDUCER 


Features  &  Applications 

•  Excellent  long-term  stability! 

•  Fully  compensated  and  signal 
conditioned! 

•  Optional  cable  available  for  low-cost 
stand  alone  applications. 

•  Great  for  applications  such  as  fume 
hoods,  clean  rooms,  VAV,  duct  static, 
medical  applications  and  many  more! 

•  Can  be  used  on  some  water  applications! 

•  Virtually  no  position  sensitivity! 

•  Custom  ranges  available! 

•  Standard  2-year  warranty! 

•  Direct  to  board  solder  connection 
available  for  easy  integration  into  most 
OEM  products. 

For  the  Most  Recent  Product 
Information  Visit  Our  Website  at 
www.autotraninc.com 


Specifications 

Accuracy* 

Stability 

Thermal  Effects:  (zero) 
0.5”  Range 

Thermal  l- fleets  (span) 

Overpressure: 

Pressure  Range 

Compensated  Range. 
Media: 


Operating  Humidity: 
Operating  Temp 
Input  Supply : 

Supply  Current 
Load  Resistance: 
Output  Signal. 


Adjustments: 
electrical  Connections: 

Pressure  Connections: 
Dimensions: 


+/-  1%  F.S.O. 

+/-  0.5%  F.S.O./yr. 

+/-o.o75%  f.s.o  .rc 
(+/-0.042%F.S.O./°F) 

+/-0. 1 50%  F.S.O./°C 
(+/-0.083%F.S.O./°F) 

+/-0.005%  F.SO/°C 
(+/-0.003%F.S.O./°F) 

20  PS1  or  2x  FSP  whichever  is 
greater 

0.5”  WC  to  40  PSI  (Custom  & 
Bi-directional  ranges  available) 

10°  to  5CPC  (50°  to  122°F) 

Limited  only  to  media  that  will  not 
attack  Polypheny  lene  Sulfide  (PPS). 
Polyetherimide  (PEI).  Silicon,  or 
Fluorosilicone,  Silicone  RTV.  Note 
that  liquids  are  allowed  in  either  or 
both  ports. 

90%  RH  non-condensing 
-25°C  to  70°C  (-13  to  I58°F) 

12-24  VDC 
<  10mA 
2KQ  minimum 

1-5  VDC,  I -6 VDC.  1-10VDC** 
(call  for  custom  voltage  outputs, 
most  available) 

Offset  60%  of  F.S.O.  minimum 
0.025”  square  pins  on  0. 1  ”  centers 
(5)- Board  Mountable*** 

Barbed  fitting  for  1/8”  I  D.  tubing 
Approx.  1 .50”  x  1.35”  x  0.9” 
(3.81cm  x  3.43cm  x  2.29cm) 


*  Includes  non-linearity,  hysteresis,  and  non-repeatability  at  a  fixed 
temperature 

**  I8-24VDC  input  required  for  1-10VDC  output 
***  Direct-to-board  solder  connection  version  available  upon  request, 
also  optional  cable  available 

AutoTran 

INCORPORATED 

7424  Washington  Avenue  South  •  Eden  Prairie.  MN  55344 
(800)  735-8998 •  (952)  942-8743 •  FAX  (952)  942-8753 
www.autotraninc.com 


ID  Number  1210-137 


Page  I  of 2 


Revision  C 


Issue  Date:  2/26/02 
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How  to  Order 


A  typical  order  number  consists  of  the  model  number,  type,  pressure  range,  input  supply,  and  output  signal. 


Model  T> ne  Pressure  Range  Input  Supply 

850  D  Differential  User  desired  range  12D  12-24  VDC 

or  Gage*  from  0/0.5”  WC  -  0/40  PSI  24D  =  1 8-24  VDC 

(Bi-directional  ranges  available) 


Output  Signal 

0  1-5  VDC 

1  =  1-6  VDC 
2=  1-10  VDC** 

5  =  custom  output 


Example:  850D  5”  WC  12D4  (Model  number  850.  differential  or  gage  type  unit,  pressure  range  of  0-5” 

WC.  input  supply  of  12-24  VDC,  output  signal  of  1-5 VDC) 

*  With  a  gage  application,  one  of  the  two  pressure  ports  is  vented  to  atmosphere. 

**  Requires  18-24  VDC  input 

Optional  Items 
ABS  plastic  enclosure 

Warranty 

Any  regular  Auto  Tran  product  that  proves  to  be  defective  in  material  or  workmanship  within  a  warranty  period 
of  two  (2)  years  will  be  repaired  or  replaced  when  returned  freight  prepaid. 

.Note:  Auto  Tran  reserves  the  right  to  change  its  specifications  at  any  time,  without  notice.  Auto  Tran  is  not 
expert  in  the  customer’s  technical  field  and  therefore  does  not  warrant  the  suitability  of  its  product  for  the 
application  selected  by  the  customer. 
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